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Abstract. A convergence theorem is proved for a class of Nystrom methods 
for weakly singular integral equations on surfaces in three dimensions. Fred- 
holm equations of the second kind as arise in connection with linear elliptic 
boundary value problems for scalar and vector fields are considered. In con- 
trast to methods based on product integration, coordinate transformations and 
singularity subtraction, the family of Nystrom methods studied here is based 
on a local polynomial correction determined by an auxiliary system of moment 
equations. The polynomial correction is shown to remove the weak singular- 
ity in the integral equation and provide control over the approximation error. 
Convergence results for the family of methods are established under minimal 
regularity assumptions consistent with classic potential theory. Rates of con- 
vergence are shown to depend on the regularity of the problem, the degree of 
the polynomial correction, and the order of the quadrature rule employed in 
the discretization. As a corollary, a simple method based on singularity sub- 
traction which has been employed by many authors is shown to be convergent. 



1. Introduction 

The method of integral equations has a long and rich history in both the analysis 
and numerical treatment of boundary value problems [7J |H [T3J [T51 [55J 1231 [13 HH 
[3H1HS] ■ On the analysis side, the method provides a classic approach to the study of 
existence and uniqueness questions. On the numerical side, the method provides an 
alternative formulation which can be efficiently discretized. In the basic approach, a 
boundary value problem described by a partial differential equation on a domain of 
interest is reduced to an integral equation on the bounding surface. The unknown 
scalar or vector field throughout the domain is represented in terms of one or more 
integral potentials which depend on an unknown surface density. The representation 
can take various different direct and indirect forms and the integral potentials 
usually involve kernels that are at least weakly singular. The reduction in dimension 
makes the method of integral equations extremely attractive, especially for problems 
on exterior domains. Indeed, the method has been used in many different modern 
applications in acoustics and electromagnetics [21 [33], hydrodynamics [HJ 03], 
elastostatics [U [37], and molecular modeling [Tl [51 [51 flu] [50] . 

Once a problem has been reduced to an integral equation, various different types 
of methods are available for its numerical resolution [7JIT8J123, 30 , 45j. In Galerkin 
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methods, the unknown density and integral equation are projected onto a finite- 
dimensional subspace of functions on the surface, where the projection is defined 
through an inner-product on the space. Basis functions for the space may be piece- 
wise polynomials, splines or wavelets with local supports [H [IH [34], or related 
eigenfunctions with global supports [6] HU [TT1 [21] . A method is completed through 
the specification of appropriate quadrature formulae for the evaluation of double 
surface integrals in the inner-product. In collocation methods, the unknown density 
and integral equation are projected similar to before, but the projection is defined 
through pointwise interpolation on the surface. As a consequence, only single sur- 
face integrals arise, which is computationally attractive. In Nystrom methods, the 
integrals in the integral equation are directly approximated by a quadrature rule, 
which leads to a direct approximation of the unknown density at the quadrature 
nodes, which can then be interpolated over the surface. For problems with contin- 
uous kernels, the Nystrom approach leads to simple, efficient and well-understood 
schemes, and for problems with weakly singular kernels, various modifications are 
required to handle the singularities, such as product integration, coordinate trans- 
formations and singularity subtraction [4| 17] [TT | fT2 | [T8 | [30 t l36 l H4 l [47 ] . 

In this article, we study a class of modified Nystrom methods based on local 
polynomial corrections. The methods studied here share some similarities with 
product integration and singularity subtraction methods, but have various impor- 
tant differences. The basic approach is to employ standard quadrature weights 
away from the diagonal of the weakly singular kernel, and generalized or corrected 
weights near the diagonal. In product integration methods, the corrected weights 
are defined so as to include the action of the kernel on a mesh-dependent basis of 
interpolating polynomials. In the methods considered here, the corrected weights 
include the action of the kernel on a mesh-independent basis of local floating poly- 
nomials, and furthermore the weights themselves are expressed as the pointwise 
values of such a polynomial. The local polynomial that generates the corrected 
weights is determined through an auxiliary system of moment equations. Such a 
local polynomial correction is shown to implicitly regularize the weak singularity in 
the kernel, as would occur in a singularity subtraction method, and provide control 
over the approximation error. 

We establish a convergence theorem for a family of modified Nystrom methods 
for integral equations on surfaces in three dimensions. We restrict attention to 
Fredholm equations of the second kind with weakly singular kernels as arise in con- 
nection with linear elliptic boundary value problems for scalar and vector fields. In 
contrast to methods based on product integration, coordinate transformations and 
singularity subtraction, we consider methods based on a local polynomial correc- 
tion as described above. Convergence results in the standard maximum norm are 
established under minimal regularity assumptions consistent with classic potential 
theory. Rates of convergence are shown to depend on the regularity of the surface 
and data, the degree of the local polynomial correction, and the order of the un- 
derlying quadrature rule. In the minimal regularity case, there is no lower bound 
on the rate, and in the smooth case, there is no upper bound on the rate. Indeed, 
arbitrarily high rates of convergence can be achieved with local polynomial cor- 
rections of arbitrarily high degree. Various important assumptions on the weakly 
singular kernel and the quadrature rule are discussed in detail. As a corollary of 
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our result, we show that a simple method based on singularity subtraction which 
has been employed by many authors is convergent under the same assumptions. 

The mathematical theory of convergence of Nystrom methods of integral equa- 
tions of the second kind on regular surfaces is a well-studied subject. The conver- 
gence theory for standard methods on problems with continuous kernels is classic 
I18[ 130] . Similarly, the theory for product integration methods on problems with 
weakly singular kernels is also well-established [3 [HI [30] , although some of the 
hypotheses may be difficult to verify in three-dimensional problems. Convergence 
results for methods based on floating polar coordinate transformations are described 
in [11] [31] [47] • Convergence results for methods based on the classic idea of singular- 
ity subtraction as introduced by Kantorovich and Krylov [27] have been established 
for one-dimensional problems in [4], and have been studied for higher-dimensional 
problems in [44j . Methods similar to those introduced here were considered in 
[12] . but no mathematical convergence results were given. Numerical experiments 
which illustrate several aspects of the lowest-order method of the family introduced 
here can be found in [36j . Related convergence results under different regularity 
assumptions are given in |35) . 

The presentation is structured as follows. In Section 2 we define the class of 
boundary integral equations to be studied. We outline various important assump- 
tions on the surface and kernels associated with the equation, and introduce no- 
tation and results that will be needed throughout our developments. In Section 3 
we define our family of Nystrom methods with local polynomial corrections as de- 
scribed above. We outline essential assumptions on the quadrature rule and other 
elements of the method and then state our main result. In Section 4 we provide 
a proof of our result. The proof is based on the theory of collectively compact 
operators and several technical lemmas; the latter are established herein and may 
be of independent interest. 

2. Boundary integral equation 

2.1. Preliminaries, notation. Throughout our developments, we consider a sur- 
face F in M 3 consisting of a finite number of disjoint, closed, bounded and orientable 
components. We invoke the standard assumptions of classic potential theory |2"2"ll2"8"] 
and assume that r is a Lyapunov surface. Thus: 

(LI) there exists a well-defined outward unit normal v(x) and tangent plane T X T 
at every x £ T, 

(L2) there exists constants C > and < A < 1 such that Q(v(x),v(y)) < 

C\x — y\ x for all 1,1/ £ F, where 9(y(x),v(y)) is the angle between v{x) 

and v(y), and \x — y\ denotes the Euclidean distance, 
(L3) there exists a constant d > such that, for every x G -T, the subset r n 

B(x, d) is a graph over T X T, where B(x, d) denotes the closed ball of radius 

d centered at x. 

We refer to A and d as a Lyapunov exponent and radius associated with r, and 
for any x £ r, refer to B(x, d) as a Lyapunov ball at x. Notice that, if (L2) and 
(L3) hold for some values of A and d, then they also hold for all smaller values. 
We assume that values for A and d are fixed once and for all, and for simplicity we 
assume from the outset that A = 1. 

For any xq £ r, we use r x0i d to denote the portion of r within the Lyapunov ball 
at xq, and use f2 Xo ,d to denote the image of r xo< d on T Xo r under projection parallel 
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to is(xo). We refer to r xo ,d as the Lyapunov patch at xq. Without loss of generality, 
we identify fi Xo .d with a subset of R 2 , and identify xq with the origin. We reserve 
the notation T Xo T to indicate the tangent plane considered as a subspace of W 3 . 
The Lyapunov condition (L3) implies that the map 

(1) y = i>x (Q, y&r XOid) £en xo . d , 

defined by projection parallel to v(xq), is a bijection. We refer to y = ipx (0 with 
inverse £ = 4> x ^{y) as a local Cartesian parameterization of r at xq. When there 
is no cause for confusion, we abbreviate ipxoiO by y Xo (t;), and abbreviate if>~ }(y) 
by £x (z/)- For any xq € T and y G r x0td , the local Cartesian coordinates £ Xo (y) 
are uniquely defined up to the choice of orthonormal basis {£x ,i> ^0,2} in T Xo r. 
Specifically, we have 

( 2 ) 6* = ^xo,aiv) = iV-Xo) -t XOia , a = 1,2. 

Throughout our analysis, we will also have need to consider a local polar param- 
eterization of r at xq of the form 

(3) i, = c lar (p,6 = i M^1), ycr X0 , d , Mefiffi, 

where ^° d T is a subset of R + x S. Here R+ denotes the non-negative reals and 
S denotes the unit circle. Notice that the map £ = m(p,£) = pt; is one-to-one for 
all £ ^ 0, with Jnverse (p, £) = n7 _1 (£) = (|£|,£/|£|), and is many-to-one at £ = 0, 
with = w(0, £) for all ( g S 1 . For any Xq (z T and 7/ G -T^o.d with y ^= xq the map 
in ([3]) is invertible. Specifically, for any choice of orthonormal basis {t Xo ,i,t Xo ,2} m 
T Xo r, we have 

(4) (p,£) = [^° la r 1 (y) where p=\y-x \, t = ( f~ ^ , a = 1, 2. 

We use C m (U, V) to denote the space of m-times continuously differentiable 
functions from U C R" into V C M fe , and C m,1 ([/,F) to denote the subspace of 
functions whose derivatives up to order m are Lipschitz. When U is not an open 
set or the closure of an open set in R™, we interpret C m (U, V) to mean functions 
possessing an m-times continuously differentiable extension to an open set or the 
closure of an open set which contains U . For instance, when U is a Lyapunov patch 
on r, we may consider extensions that are constant in the normal direction. We say 
that r is of class C m , and use the notation r G C m , if ip xo G C m ((2 Xo . d ,R 3 ) and 
Tp^ 1 G C m (r x0id , R 2 ) for every x G r. Similarly, we say that r is of class C™' 1 , and 
use the notation T G C* m '\ iftp xo G C m ^{Q Xo , d , R 3 ) and ^ G C" 1 - 1 (T^d, R 2 ) for 
every £0 G -T, and additionally the Lipschitz constant for each derivative is uniform 
in xo. 

To any / G C 1 ^, R fe ) we associate a surface derivative Df G C°(T, R fex3 ). 
Specifically, for any Lyapunov patch r X(ud and any point y = ipx (0 G r x0jd , we 
define Df(y) by 

d( f°ip Xo ) w _ 

(6) £/(y)^ = o, V^G^r^. 

Since any vector v can be decomposed as v — Vt + v v , where V( e T y r and 
v v G [T,,/ 1 ^, we have Df(y)v = Df{y)v t . Thus L>/(y) G R kx3 is defined for all 
vectors v G R 3 . For arbitrary / G C 1 (_T, R fc ), the matrix Df(y) is equivalent to 



(5) D}{y)t= — p, Vt = —^p£T y r, pE 
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the usual derivative matrix of an extension of /, where the extension is constant 
in the normal direction to r. When / £ C 1 (r,M. k ) is the restriction of some given 
g £ C 1 (R 3 ,K' C ), the matrix Df(y) is equivalent to the composition of the usual 
derivative matrix of g with the projection matrix for orthogonal projection onto 
T y T '. In the case when g is constant in the normal direction, the projection matrix 
can be replaced by the identity matrix and the two descriptions coincide. Notice 
that, for any curve y(r) = V ; a;o(?( r )) m an y Lyapunov patch r xo ,d, we have 

(7) ^f(y(T))=Df(y(T))±y(T). 

Given a function of the form f(x, y) with x, y £ r, we use the notation f x (y) and 
f v (x) to denote the functions obtained by fixing the values of x and y, respectively. 
Similarly, we use the notation D x f(x, y) and D y f(x, y) to denote the surface deriva- 
tives of f{x, y) with respect to x and y, respectively, where the other variable is held 
fixed. Hence D x f(x,y) = Df y (x) and D y f(x,y) — Df x (y). Higher-order surface 
derivatives for a function f(y) are defined in a natural way. For instance, if we iden- 
tify R kx3 with R 3k and Df £ C 1 (r,M 3k ) 1 then D 2 f = D(Df) £ C°(T, R 3A;x3 ) is 
defined exactly as above. We denote the diagonal subset of r x r by T = {(x, y) £ 
r x r | y = x). Moreover, for any 6 > we also define an open 6- neighborhood of 
the diagonal by Ys = {(x, y) € r x r \ \y — x\ < 5}. Throughout our developments, 
we use | ■ | to denote a Euclidean norm or the measure of a surface, as determined 
by the context. 

2.2. Problem statement. Given a Lyapunov surface r in three-dimensional space 
R 3 , we consider the problem of finding a function tp £ C°(r,M. k ) that satisfies 

(8) ap-Aip = f, 

where c ^ is a given constant, / £ C°(r,R k ) is a given function, and A : 
C°(T, R k ) —5- C°(r, R fe ) is a given integral operator of the form 

(9) Atp^S^ + Kip, 
where 

(10) (9<p)(x) = J r G(x, yMy)dA y , pi<p)(x) - J H(x, yMy)dA y . 

In the above, dA y denotes an infinitesimal area element at y, G £ C°(r x F, M' cx ' c ) 
is a continuous kernel, and H £ C°(r x r\T, M. kxk ) is a weakly singular kernel, 
which is unbounded along the diagonal subset T. The above equation is typically 
considered with k — 1 or 3. The first case corresponds to a scalar- valued problem 
in three-dimensional space, as would arise in applications to electrostatics, and the 
second case corresponds to a vector-valued problem in three-dimensional space, as 
would arise in applications to elastostatics and hydrodynamics. 

Our study of ([8]) will rely on various assumptions about the surface r and the 
kernels G and H. Specifically, we assume: 

(AO) r £ C m+1,1 for some m > 0, with Lyapunov radius d > 0. 
(Al) G £ C m ^(r x r,R fexfe ). 

(A2) H £ C m ' l {r x r\T s ,R kxk ) for every < S < d. Moreover, H can be 
decomposed as 

(11) H(x,y) = u(x,y)/\x-y\ 2 -^ 
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for some function a e x r\Y,R kxk ) and exponent < /x < 1. Fur- 

thermore, on r x r\T the component functions u%j (1 < i,j < k) satisfy, 
for some positive constant C, 

(12) \ Uij (x,y)\<C, 

(13) \D s x Uij(x,y)\ < C\y-x\~ s , l<s<m+l, 

(14) \D a yUij (x,y)\ < C\y-x\- s , 1 < s < m. 

(A3) When m > 1 we additionally assume that /i = 1 and that u has the prop- 
erties 

(i) < olar eC ' 1 (i2P ol f,M fcxfe ), Vx a er, 

(h) < lar (0,-a = < lar (0,a VCeS, 
(hi) ^eC^ 1 ^,,!'^), Vx er, 

where Lipschitz constants are uniform in xo and u£° lar = u Xo o ip^.° lar and 

(15) u A (y) = {[< lar ^o-< ar (o^li ( M>^r]-u,)' 

[0, y = x . 

Assumption (AO) states that r must be at least of class C ' , which implies 
that -T is differentiable, with a curvature that is defined almost everywhere and 
bounded. Assumption (Al) states that the regularity of G may be one order less 
than that of r, which will be convenient for our analysis. Assumption (A2) is 
essentially the definition of a weakly singular kernel. Typically, the parameter /i is 
identified with the Lyapunov exponent A, but it is not necessary to do so. Away 
from the diagonal, the regularity of H is assumed to be the same as that of G. 
Assumption (A3) outlines additional regularity conditions on H in the case when 
r is at least of class C 2 ' 1 . The most specific condition is (A3)(ii), which states that 
the function w!j!° lar is an even function on S at p = 0. As we will see, this implies 
that certain local moments of H will vanish. This condition will be important in our 
analysis of numerical methods and is related to the classic Tricomi condition that 
arises in the study of singular integral operators [46 a . This condition is satisfied 
by weakly singular kernels H arising in different applications, for example, the 
classic double-layer kernel in three-dimensional electrostatics, and incompressible 
elastostatics and hydrodynamics. We remark that condition (A3) is required for 
only the lowest-order method studied herein corresponding to a local polynomial 
correction of degree zero; it is not required for the higher-order methods. 

2.3. Solvability theorem. The following result establishes the solvability of the 
boundary integral equation in . Its proof follows from the classic Frcdholm The- 
orems for compact operators [30l [38] and is omitted for brevity. The compactness 
of A = S + 3~C will be established in Lemma 14.51 



Theorem 2.1. Under conditions (A0)-(A2) the operator A on C°(r, R fe ) is com- 
pact. Hence, provided c is not an associated eigenvalue, there exists a unique so- 
lution ip e C°(r,R k ) for any closed, bounded Lyapunov surface r G C ' and 
boundary data f £ C°(r,R k ). 

Thus, under mild conditions, the boundary integral equation in ([8]) has a unique 
solution if for any data /. Moreover, because the operator cl — A can be shown to 
have a bounded inverse, it follows that ip depends continuously on /. The regularity 
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of the solution ip is determined by that of the data / and properties of the operator 
A. When A is a smoothing operator, so that Aip is smoother than (p, it follows from 
that if is in the same regularity class as /. The form of the operator A arises in 
various applications as mentioned above. The component "K may be interpreted as 
a double-layer potential. The component S may be zero, or may be interpreted as a 
range completion term necessary for the uniqueness of solutions. In the latter case, 
S may represent a potential due to a discrete or continuous distribution of sources 
located away from the surface. Such operators can be motivated beginning from a 
generalized boundary integral representation of a field by a linear combination of 
single- and double-layer potentials, and then moving the single-layer potential to 
an offset surface, and then possibly also shrinking the offset surface to a curve or 
point, or collection thereof. Various forms of such operators have been previously 
considered [19l[22l[2ll[29l|33l|38l|40l|4T]. 



3. NYSTROM APPROXIMATION 

3.1. Mesh, quadrature rule. We consider an arbitrary decomposition of r into 
non-overlapping quadrature elements r e , e = l,...,E, each with area \r e \ > 0. 
To any such decomposition we associate a size parameter h = max e (diam(-T e )) > 
0. For simplicity, we assume that the elements are either all quadrilateral or all 
triangular. In our analysis, we consider sequences of decompositions with increasing 
E, or equivalently, decreasing h. We will only consider sequences that satisfy a 
uniform refinement condition in the sense that the area of all elements is reduced 
at the same rate. Specifically, we assume 

Ch 2 < \r e \ < C'h 2 , Ve = l,...,£, E>E . (A4) 

Here C, C and Eq are positive constants whose values may change from one ap- 
pearance to the next. 

In each element r e , we introduce quadrature nodes x e q and weights > 0, 
q = 1, . . . , Q, such that 

. E r e q 

(i6) / f{ X ) dA x = Y, /(*) ^*EE /K)^ e - 

Jr e=l Jre e=l 9=1 

Without loss of generality, we assume that the quadrature nodes and weights are 
defined by mapping each element r e to a standard, planar domain using a local 
parameterization of the same regularity as the surface, and applying a local quadra- 
ture rule in the standard domain. In this case, the Jacobian of the parameterization 
would be included in the weights W q . We assume that the quadrature weights re- 
main bounded and that the quadrature points remain distinct and in the element 
interiors. Specifically, for any sequence of decompositions that satisfy the uniform 
refinement condition, we assume 

Ch 2 < E 9 Q =i% e < C'h 2 , Ve =!,...,£, E>E , 
Ch<dist(x e q ,dr e ) < h, Vq = l,...,Q, e = l,...,E, E > E 0l 
Ch< min \x e -x e ',\<h, Vq,q' = 1,...,Q, e, e' = 1, . . . , E, E > E . 

(e,q)#(e',q') 

(A5) 
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To quantify the error in a quadrature rule for a given function / on a given 
surface r, we introduce the normalized local truncation errors 



(17) T ( e ,/,fc) = _L 



r 9 

/ f(x)dA x -J2f(x e q )W t 
J r* q=1 



l,...,E. 



For sequences of decompositions that satisfy a uniform refinement condition, we 
require that the above truncation errors vanish uniformly in e depending on prop- 
erties of /. Specifically, we assume there exists an integer £ > 1 such that 

r(e, /, h) — > as h — > uniformly in e, for each / £ C°(r, M), 

r(e, /, h) < Ch l uniformly in e, for each / e C l - l ' l {r, R). ' A<> ' 

In the above, the constant C may depend on /, but is independent of e and h, and 
the integer I > 1 is called the order of convergence of the quadrature rule. 

For convenience, we will often replace the element and node indices e = 1, . . . , E 
and q = 1, . . . , Q with a single, general index a = 1, . . . , n, where n = EQ. We will 
use the multi- and single-index notation interchangeably with the understanding 
that there is a bijective map between the two. 

3.2. Partition of unity functions. To each quadrature node x a in a decompo- 
sition of r we associate nodal partition of unity functions ( a ,Ca £ C°(-T, R). We 
assume that these functions take values in the unit interval and are complementary 
so that their sum is equal to one. Moreover, we assume that £ a vanishes at least 
quadratically in a neighborhood of x a , and that the support of Q a is bounded from 
above by a multiple of the mesh parameter h. Specifically, we assume 

< C«(aO> Cafc) < 1, Ca(x) + ( a (x) = 1, Vxe r, 

C\x - x a \ 2 - (^ 7 ) 

\(a(x)\ < ^ 2 a , diam(supp(Ca)) < Ch, Va = 1, . . . , n, n > n . 

The functions £ a and £ a will play an important role in the family of numerical meth- 
ods that we introduce and in the associated convergence proof. Specifically, these 
functions will help isolate the weak singularity of the kernel H at each quadrature 
point and control the numerical error there. 

3.3. Discretization of integral equation. Let a decomposition, quadrature rule, 
and nodal partition of unity functions for r be given. For any function /, let S/t 
and yih denote approximations to S and !K of the form 

n n 

(18) (5 h .f)(x) = Gb(x)f(x b ), (X h f)(x) = J2 H b (x)f(x b ), 

h=l b=l 

where and H b are functions to be specified. Moreover, in view of ©, let ifih 
denote an approximation to tp defined by 

(19) aph - A h cp h = f where A h = 5 h + Ji h . 

The equation for ip^ can be reduced to an n x n linear system for the nodal 
values (fih(x a ). Indeed, from (TT51) and we get, for each a = 1, . . . , n, 

n n 

(20) c(p h (x a ) - ^ G b (x a )ip h (x b ) - ^ H b (x a )(p h (x b ) = f(x a ). 

6=1 b=l 
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Notice that, for every solution of the discrete system (|20|) . we obtain a solution of 
the continuous system (|19|) . namely 

n n 
b=l 6=1 

Moreover, the converse is also true; every solution of the continuous system provides 
a solution of the discrete system by restriction to the nodes. 

The method is completed by specifying the functions G b (x) and H b (x). In view 
of (|TU|) . (flo) and (fT5j) . we define Gb{x) = G(x,x b )W b . However, due to the weakly 
singular nature of the kernel function H(x,y), a similar definition cannot be made 
for Hb{x), because H (x, x b ) is not continuous or even defined when x — x b . Instead, 
for any integer p > we define 

(22) H b {x) = ( b {x)H{x,x b )W b + ( b {x)R x (x b ), 

where R x {x b ) is a local polynomial of degree p at x, evaluated at quadrature point 
x b . By a local polynomial at x we mean a polynomial in any system of rectangular 
coordinates in the tangent plane with origin at x\ see (|24p below. For a Lyapunov 
surface r, such polynomials are well-defined in a neighborhood of each point x, and 
there is a uniform bound on the size of this neighborhood; indeed, they are well- 
defined in the Lyapunov patch r x c i- The unknown coefficients in R x are defined 
by enforcing the local moment conditions 

(23) (H h rj x f x )(x) = (Jtr) x f x )(x), for all local polynomials f x up to degree p. 

Here rj x € C m,1 (r, [0, 1]) is any given cutoff function which is identically one in a 
fixed neighborhood of x, for example r x ^u, and which is identically zero outside 
some fixed neighborhood of x, for example r x d / 2 - The local polynomials R x and f x 
can be described as floating since they are defined in a tangent plane that depends 
on x. 

The basic structural form of H b is similar to G b , but with a correction in a 
neighborhood of x b . By construction, for all x outside a neighborhood of x b , we 
have Cb — 1 and ( b = 0, so that H b is defined by the quadrature rule applied 
to H. On the other hand, as x approaches x b , we have Cb - ► and Cb 1, so 
that H b is determined by the local polynomial R x . The value of R x at a node x b 
can be interpreted as a generalized quadrature weight as would arise in a product 
integration method. There are various freedoms in the choice of r] x used in the 
moment conditions. For the lowest-order method with p = 0, the local polynomial 
R x reduces to a constant polynomial which can be extended to the entire surface, 
and the cut-off function is unnecessary and can be taken as unity. However, for 
higher-order methods with p > 1, the local polynomial R x in general cannot be 
extended to the entire surface, and a non-trivial cut-off function as described above 
is necessary. Moreover, a different cut-off function could be used for each different 
coefficient in R x . For convenience, we employ a single cut-off function r\ x for all the 
moment conditions in our analysis below. 

The moment conditions in (1231) lead to a linear system of equations for the 
unknown coefficients of the local polynomial R x , which can be solved for any given 
point x. Specifically, for any given leJ 1 , the local polynomial R x : r X: d K fexfe 



(21) <p h {x) = - 
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of degree p > has the form 
(24) 

where the usual summation convention on pairs of repeated indices a\, . . . , a p is 
implied. Substituting (j2"4"]l into (l2"2"j) . and using the standard basis for the local 
polynomials f x , we find that the moment conditions in (|23p lead to a linear system 
for the k x k coefficients C x ,o & n d Ca;. Ql through Cx, ai a 2 ---a p , namely 

p 

M^fiCxfi + M x,a 1 ---a s Cx,a 1 ---a s = A°, 
s=l 
P 



(25) 



s=l 

When p = 0, the indicated sums and all equations except the first are empty and 
the system reduces to the single equation M° C Xl o = A°, where 



Ml Q = ^2Cb{x)vM GR, 



(26) ' 1 



A°= / ff(x,tf)»fe(tf)di4 J ,-VC6(a:)fr(a!,a!6)tb(aj6)W6G 



6=1 

When p > 1, all sums and equations are non-empty and the system has the full 
form indicated in (|25|) . In addition to M® and A°, we have, for 1 < s, t < p, 

n 

M x, ai - as =^2(b{x)Vx{Xb)tx,a ll (Xb) ■ ■ ■ Zx,a s (xb) GR, 
6=1 

n 

(27) M^ - ft = ^C&teKfo&.ftfo) • ' • £x,p t {%b) e R, 

6=1 

n 

M x]Z-<x a = ^2 Cb(x)r]x(Xb)^x,ai(Xb) ■ ' ' (,x, as {xb)£,x,f}A x b) ' ' ' £x,P t ( X b) G R, 
6=1 



and 



(28) 



H{x,y)r) x (y)£x,/3i(y) ■ ■■(x,p t (y) dA y 



G M fcxfc . 



- ^C 6 (a;)iT(£, x b )rj x (x b )^ x ^ 1 (xb) ■ ■ ■ £,x,p t (x b )W b 

6=1 

In contrast to what the number of indices would suggest, the number of indepen- 
dent k x k equations in (1251) grows only quadratically in p. This follows from the fact 
that all quantities in the equations are fully symmetric in the indices ot\ ■ ■ ■ a s and 
/3i • • • fa for all 1 < s, t < p. Accounting for symmetry, the number of independent 
components of C x , ai ---a s is not 2 s , but only (s + 1). Summing over s = 1, . . . ,p, 
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and including the single contribution from C x ,o, we find that the number of inde- 
pendent k x k unknowns, or equivalently, independent k x k equations, is equal 
to (p + l)(p + 2)/2. Notice that the implementation of the system in ([2"5"]) requires 
various integral moments of the weakly singular kernel H. These moments can be 
evaluated numerically using techniques such as Duffy |15j or local polar coordinate 
[TTj l3T[ [47] transformations. In the case when p = and the cut-off function is 
taken as unity, the required moment is known analytically in many applications. 
However, in the case when p > 1 and the cut-off function is non-trivial, the required 
moments must generally be obtained numerically. 

For any given point x G r, the solvability of the linear system in ([25]) for the 
coefficients of R x depends on the support of the nodal functions £&. Specifically, let 
Jx = {b | Cb{%) > 0} and consider sufficiently refined surface decompositions such 
that rj x {xb) = 1 for all b G J x . For any local scalar- valued polynomial f x , let F x 
denote the vector of size (p + + 2)/2 of independent coefficients, where each is 
weighted according to its multiplicity. Moreover, let M x denote the corresponding 
square, symmetric coefficient matrix of size (p+ l)(p + 2)/2 associated with the 
independent equations in ([25]) . Then by direct computation we find 

n 

(29) F x -M x F x = J2(b(x)[f x 

6=1 beJ* 

From this we can deduce sufficient conditions for the positive-definiteness of M x and 
hence the unique solvability of (f2"5j). Specifically, for each x G r, it is sufficient that 
J x be non-empty, and that the only polynomial of degree p which satisfies f x {xb) = 
for all b G J x be the zero polynomial. In view of the number of polynomial 
coefficients, this condition implies that J x must contain at least (p + l)(p + 2)/2 
quadrature nodes for each icf. When x is itself a quadrature node, this implies 
that the support of each function must contain at least as many nodes. Hence 
the support of the functions determines the solvability of ([2"5]) . In our analysis, 
we will assume that M x has a uniformly bounded inverse, namely 

\M~ l \<C VrrGT, n>n . (A8) 

A straightforward choice of functions which satisfy conditions (A7) and (A8) in 
the case p = is described below; see also [36] . 

The implementation of the numerical method is centered upon ([20]) . ([2Tj) and 
([2!))) . For each quadrature node x a , the linear system ([23|l is solved to obtain the 
local polynomials R Xa , which are then used in the linear system ([20]) to obtain the 
nodal values tph(x a ). Once these nodal values are determined, they can be extended 
to a continuous function. Specifically, given any point x, the linear system in ([25]) 
can be solved to obtain the local polynomial R x , which can then be used in the 
interpolation equation ([2"T]) to obtain the value of (fih(x). Notice that, in general, the 
construction of the local polynomials R Xa and R x requires the evaluation of local 
moments of the weakly singular kernel H, and the evaluation of local Cartesian 
coordinates ^ Xa and £ x . As will be shown later, the values of R Xa and R x are 
independent of the choice of orthonormal basis associated with £ Xa and £ a . 

3.4. Illustrative example. Here we illustrate the form of the nodal equations ([20]) 
in the lowest-order case with p = 0. In this case, the equations take a particularly 
simple form and the evaluation of local Cartesian coordinates is not necessary, and 
the evaluation of weakly singular integrals is typically not necessary. 
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We begin by describing partition of unity functions Ca and Ca which satisfy 
conditions (A7) and (A8). Consider an auxiliary decomposition of r into Voronoi 
cells, where each cell contains a single quadrature point x a . Each cell can be 
mapped to a unit circle, with x a mapped to the center. A simple quadratic function 
z = (x 2 + y 2 )/2 in the unit circle can be mapped back to the cell as the central part 
of Ca- We can then introduce an offset boundary which is displaced outward from 
the cell boundary by a distance of e/2, where e = min a7 ^ \x a — Xb\- The mapped 
quadratic function, which by design has the value 1 /2 on the cell boundary, can be 
extended to achieve a value of 1 on the offset boundary, and then further extended 
to the rest of r with the constant value 1. Notice that the functions Co and 
Ca so constructed have the convenient nodal property that ^ a (xb) = 1 — 5 a b and 
(a(xb) — Sab- Moreover, the supports of £ a overlap on the entire surface and we 
have J2 a CaO) > 1/2 for all x e T. 

For the method with p = 0, the local polynomial R x reduces to a constant 
polynomial R x {z) = C x ,o, and the cutoff function r\ x can be taken as unity. The 
linear system in (|25|) reduces to the single equation M® C Xt o = A°, which implies 

(3Q) c = J r H ( x >y) dA v ~ T,b=i(b{x)H(x,x b )W b 

x '° V^ n T ! \ 

2^b=i(b{x) 

This solution is well-defined and bounded for any sequence of decompositions by 
properties of the nodal functions Ca and Ca and the kernel function H . In var- 
ious applications, the weakly singular integral in the above expression is known 
analytically and hence numerical evaluation is not necessary. 

The nodal equations in (|2"01) take a particularly simple form. Indeed, because 
our choice of the nodal partition of unity functions Ca and Ca has the property that 
(a(xb) = 1 — S a b and Ca{x b ) = 5 a b, the equations become 



l{x a )(p h {x a ) - } G(x a , x b )ip h (x b )W b 

6=1 n 

(31) -J2 H ( X a > X b)[Ph{Xb) - <Ph(x*Wb = f{Xa), 



b=l 



where j(x a ) = cl — J r H{x a , y) dA y . This discrete system is similar to the classic 
singularity subtraction method discussed by various authors [H [27]. The factor 
[<Ph(xb) — <fh{x a )] can be interpreted as cancelling the weak singularity in H{x a , x ). 
Indeed, since the sum extends over b ^ a only, the product H(x a , Xb)[<fh{ x b) — 
<Ph{%a)\ can be interpreted as vanishing when b = a. For methods with p > 1, a 
similar but higher-order cancellation can be interpreted to occur. Once the nodal 
values of tph are determined, they can be extended to a continuous function using 
the interpolation equation in (f2"Tj) . Notice that the nodal values of iph can be 
computed without explicit knowledge of the partition of unity functions. Various 
numerical examples with this method are given in [36) . 

3.5. Solvability and convergence theorem. The following result establishes 
the solvability and convergence of the locally-corrected Nystrom method defined 
in [fT8 ]) -([2"8 |) . We consider the method with a quadrature rule of arbitrary order 
i > 1, a local polynomial correction of arbitrary degree p > 0, and a surface with 
regularity index m > 0. In view of Theorem 12.11 we suppose that the constant 
c / is not an associated eigenvalue of the operator A, so that the given boundary 
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integral equation has a unique solution ip. Below we use C v to denote a constant 
depending on ip. 

Theorem 3.1. Under conditions (A0)-(A8) 7 there exists a unique approximation 
cph £ C°(r, R fe ) for any closed, bounded Lyapunov surface r € C 1 * 1 and boundary 
data f S C°(r,R fe ) for all h > sufficiently small. Moreover, if ip e C mA {r,R k ) 
and r £ C m+1 ' x , then as h 



(i) 


\\<Ph 


-<p\\ 


-+0, 


W > l,p > 0,m 


>o, 


(") 






< C v h, 


W > l,p = 0,m 


> 1, 


(iii) 


\\<Ph 






w > l,p > l,m 


> 1. 



Thus, under suitable assumptions, the method defined by (fl8 |) -(f28 | is convergent 
in the usual maximum or C°-norm. The rate of convergence depends on the order 
t of the quadrature rule, the degree p of the local polynomial correction, and the 
the regularity index m of the exact solution ip and the surface r. In the minimal 
regularity case with m = 0, there is no lower bound on the rate, and in the higher 
regularity case with m > 1, the rate is at least linear. The rate of convergence 
is independent of the order of the quadrature rule in the case when p = 0, which 
corresponds to the lowest degree of correction. Higher rates of convergence are 
obtained when p > 1, which corresponds to higher degrees of correction. Notice 
that the rates of convergence stated above are lower bounds; they could possibly be 
higher in certain circumstances, for example in smooth problems with periodicity, 
for which some quadrature rules are known to have special properties [TTJ [32] . 
We remark that the method considered here is based on open quadrature rules as 
required by condition (A5). 

As a special case, Theorem 13.11 with p = establishes the convergence of a 
method similar to the classic singularity subtraction method considered previously 
by various authors [H [27] . The two methods lead to apparently identical discrete 
systems for the nodal approximations, but differ in how the nodal approximations 
are interpolated over the surface. Convergence results for this method appear to 
be not well-known. The results derived here make crucial use of the structure of 
the nodal functions Hb(x) defined in (|22|) . the moment conditions defined in ([23)) . 
and various properties of the nodal partition of unity functions and Such 
ingredients appear to have not been considered in previous studies of the classic 
method. We remark that the linear convergence result for p — is delicate and 
relies on a Tricomi-like property of the potential !H implied by condition (A3)(ii). 
The results for p > 1 are qualitatively different and rely mainly on the regularizing 
effect of the local polynomial correction. Specifically, the linear convergence result 
in part (ii) requires condition (A3) whereas the convergence results in parts (i) and 
(iii) do not. 

4. Proof 

In this section we provide a proof of Theorem 13.11 We use the same notation 
and conventions as in previous sections. Specifically, we use C, C , C" and so on, 
to denote generic positive constants whose value may change from one appearance 
to the next, use | ■ | to denote a Euclidean norm or the measure of a surface, as 
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determined by the context, and use || • || to denote the usual maximum norm on 

c°{r,R k ). 

4.1. Collective compactness. We first outline a result, based on the theory of 
collectively compact operators [3], which plays a fundamental role in the analysis 
of Nystrdm methods [TJ [TSJ [3D] . Let A be a linear, compact operator on C°(P, M fc ) 
as given in @, and let A n (n > tiq) be a sequence of linear, finite-rank operators 
on C° (r, R fe ) as given in (|19l) , where for convenience we consider subscripts n — > oo 
in place of h n — » 0. Consider the following conditions: 

(CI) For each v G C°(r,R k ), (A n v)(x) -> (Av){x) uniformly in x G Z\ 

(C2) |(yi„w)(x)| < C for all a: G T and w G C°(T,R fe ) with ||u|| < 1. 

(C3) For every e > there exists 6 > and N > n such that \(A n v)(x) — 

(A n v)(y)\ < e for &Hx,y G T with < 5, v G C°(r,R fc ) with ||v|| < 1, 

and n > N . 

Condition (CI) states that A n converges to A pointwise in C°(r, R fc ). Conditions 
(C2) and (C3) state, respectively, that S = {A n v \ n > no, \\v\\ < 1} is an 
equibounded and equicontinuous subset of C°(r,R k ). These two conditions imply, 
by the Arzela-Ascoli Theorem, that S is relatively compact; by definition, the 
sequence A n is then called collectively compact. For such sequences the following 
well-known result holds [3 [TBI EQ] . 

Theorem 4.1. Let A and A n (n > rig) be linear operators on C°(r, R fe ), where A 
is compact and A n are finite-rank, and assume that (C1)-(C3) hold. If ' cLp — Atp = f 
is uniquely solvable for if, then there exist constants C v > and N v > uq such that 
CLp n — A n Lp n — f is uniquely solvable for ip n , and moreover 

(32) \\<p n -ip\\ < C v \\A n ip - Aifi\\, Vn>N v . 

In what follows, we establish conditions (C1)-(C3) for the operator A n = S n +!K n 
defined in (fT9|) . We concentrate on J£„ since the result for 9n is straightforward by 
continuity of its kernel. Once (C1)-(C3) are established, the result in Theorem 13. II 
will follow from a bound on ||j4 rl <y3 — A<p\\ under regularity assumptions on r and 
ip. 

4.2. Lemmata. We begin with a collection of useful results regarding the weakly 
singular kernel H . Parts (i) and (ii) below show that the integral operator defined 
by H indeed maps C°(r, M fe ) into itself, and is equibounded and equicontinuous on 
bounded subsets, hence compact. Part (iii) is an important inequality that will be 
used in the sequel. 

Lemma 4.2. Let r satisfy (AO) with Lyapunov radius d > and let H satisfy 
(A2) with exponent < fj, < 1. Then: 

(i) J r \H(x, y)\ dA y < C for all x G T. 

(ii) For every e > there exists 5 > such that J r \H(x*, y) — H(xo, y) \ dA y < e 
for all xo,x« G -T with \xq — < S. 

(iii) \H(x*,y) - H(x ,y)\ < C\x* - x \/\y - x \ 3 ^^ for all a; ,a;*,j/ G T with 
< \x„ — xq\ < d/3 and 2|x* — xo| < \y — xq\ < d. 

Proof. Parts (i) and (ii) are classic results for weakly singular integrals, see for 
example [22], where the proof of (ii) relies on (iii). For brevity, we illustrate only 
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(iii) . To begin, we consider the form of H in (fTT|) and use the triangle inequality to 
get 

u(x*,y) u(x Q ,y) 



\H(x*,y) - H(x ,y)\ 



(33) 



\y - x, 



\2-fi 



\V - s 



12— /x 



< 



\u(x*,y) - u(x ,y)\ 



\y - x. 



|2— A* 



u{x ,y) 



|y-s*| 2 |y-s | 2 ^ 



For the first term in (1331) . we consider the local Cartesian coordinate map ip Xo 
in the Lyapunov patch r x0j d, and the curve x(t) = 4>x {t£,*), < t < 1, from 
xq = V'xqCO) to = ip Xo (£*)- Using this curve, we have 



(34) 



u(x*,y) - u(x ,y) 



D x u(x(T),y)—x(T) dr. 



(35) 



\u(x*,y) - u(x ,y)\ < 



dr. 



From (fl"3|) and the relation dx/dr = (dip Xo / d£)!;* , and the facts that \dip Xo /d£\ < C 
and < C\x* — Xq\, which follow from the Lipschitz properties of ijjx and ip X g, 
we get 

f 1 C\x* - x \ 
lo \x{t) - y\ 

From the definition of x(r), we deduce that \x(t) — Xq\ is an increasing function, 
and hence |s(r) — xq\ < |x* — xq\ for all < r < 1. Moreover, by hypothesis, we 
have < |s* — s | < \y — xq\/2. These two results imply, with the help of the 
triangle inequality, 

(36) ~<^#4 re [0,1], and I < ^ < §, 

and using (f36|) in (f35j) . we find 



(37) 



|u(s*,y) - u(£ ,y)| < 



so 



\y - x*\ 

For the second term in (|33p . we notice that u| < C by fT2")) . Hence, using the 
notation r = |y — x\, we have 



(38) 



w(x ,y) 



i(xo,y) 



\y - x* 



12-Ai 



12/ - ^0 



|2-p 



< C 



1 



f 



2-p 



2-M 



c\ 



2-ui 



2-u 2-u 

r 'r* 



From ([36)) and the hypothesis < 2|x* — xo| < |y — so < d we notice that 

(39) < r„ < d and < r < 2r* < 3r < 3d. 

Regarding the factor Ir 2 .^ — r 2 ~ M | in (|38"|) , we have, since the function r 1- ' 1 is 



monotonic, 

2-u 



(40) 



r, 



o 



-2-/*| 



Jro dr 



< (2 — fx) max{r ( 



~1-U 
' ' * 



}\r* 



ro 



Combining (gOJ and {39]) with ([38]), we find 

u(x ,y) u(x ,y) 



(41) 



12/ - x, 



|2-p 



12/ - sol 



< C|r« -r | < C|x 



sol 



3-u 



2/ - s | 3 * ' 



where the last inequality follows from the definitions of ro and r* , and the straight- 
forward inequality |r* — ro| < | 

X^ Xq |. Substituting (JUJ) and (J3ZI) into (J33]), and 
again using on the first term, we obtain the desired result. □ 
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The next result can be viewed, in part, as a discrete analog of Lemma 14.21 It 
will play a central role in establishing the collective compactness of the operators 
"K n associated with the numerical method. The defining elements of the method, 
which are the surface decomposition, quadrature rule, and nodal partition of unity 
functions, are denoted by the set {r e , x^, Wg, Q, 

Lemma 4.3. Let r satisfy (AO) with Lyapunov radius d > 0, H satisfy (A2) with 
exponent < ft < 1, and {r e , x e q , Q, Q} satisfy (A4)-(A7). Then: 

(i) (,b(x)H(x, Xb) {'■= at x = Xb) is continuous in x E T for each b = 1, . . . , n 
and n > uq. 

(ii) Ylb=i \Cb(x)H(x, x b )Wb\ < C for all x G r and n > n . 

(hi) For every e > there exists S > such that J2b=i \Cb( x *)H(x*, Xb)Wb — 
(b(xo)H(x , Xb)Wb\ < e for all x Q , £ T with \x — x* \ < 5 and n > no- 
Proof. For part (i), we notice that is continuous for all x, and H{x,Xb) is 
continuous for all x ^ Xb- Hence, to establish the result, we need only verify that 
(b(x)H(x, Xb) — >• as x — > Xb- From (A2) and (A7), we get 

(42) Kb(x)H(x,x b )\ = la( " )M(x l! b)l < C \*-*»T 



x b \ 2 -f ~ h 2 

and the result follows since < fi < 1. 

For part (ii), let 5 G (0, d/3] and f3 > 10 be any given numbers and consider 
sufficiently rehned surface decompositions with E > Eg g, where Eg g > E a is a 
given integer defined such that < f3h < 5/2 for all E > Egj. We consider an 
arbitrary x = xq G r and consider collections of surface elements r e , e = 1, . . . , E, 
defined by 

(43) I Xo ,g = {e | r e C r xo , s } and I X()Sh = {e \ r e c r xthPh }, 

and note that, by design, the sets I Xo ,s, Ix ,/3h C I Xo ,S and I Xo ,g\I Xo ,gh are non- 
empty for all E > Eg t s- Moreover, we consider the decomposition 

n 

(44) J2&(xo)H(x Q ,x b )W b = § 1 (x t) )+§ 2 {x )+§3(xo), 



6=1 



where 



(45) §iM= J2 E^N)^,^)^, 



(46) s 2 (x )= E Y,Q^) H ^ x D w h 

Q 

(47) § 3 (x )= E E^^^^^r 

e^-T^o.a 9=1 

To establish the result, we show that each of the sums Si, §2 and §3 is uniformly 
bounded for all E > Egj and xq G r. Boundedness for the finite interval E G 
[E , Eg } $] is a straightforward consequence of part (i) and will be omitted for brevity. 

For the term Si in (|4"4"|) , we consider the local Cartesian coordinate map ip X() in 
the Lyapunov patch r XQt d- Since r x0t gh C r xo ,d, and ip Xo and ip X( ^ are Lipschitz, 
we have the area bound | J~Lc ,/3 | < Ch 2 . Moreover, from (A4) and the fact that 
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( U ee/x ,,jh^ e ) C r Xo ,ph, we get Ch 2 \I Xo ^ h \ < \r Xo ,ph\, where \I Xo ,ph\ denotes the 
number of elements in the index set I Xo ,0h- From these two observations, we deduce 
the uniform bound 



(48) \I X0 ,ph\<C, VE>E PtB , l0 er. 

From (|45|) we get, using (A2) and (A7), 



— , Q 

\Si(x Q )\< £ "£\e q (x )H(x ,x e q )We\, 
(4g) eeJ^u^l 

-2-^2-^ h 2 \x -x%\ 2 -» q 



eela, ,ph 9=1 

Since x q £ T e C r xo ^h we have |x — cc|| < /3ft, and from (A5) we have Y^=i W q < 
Cft 2 . Using these results together with (|4T))) and (|48p. we find that 

(50) |Si(a;o)| < Ch», VE > E 0t5 , x £ r, 

which establishes the result for §i since < fi < 1. Indeed, the sum Si(xo) —> as 
E — > oo (ft — >• 0) uniformly in xq £ r. 

For the term §2 in (|44j) . we notice first that r xo ^_i\y l C r xo ^h and that 
dist(5/ n a;0i ( (g _ 1 ) ft ,ar 2;0!/ 3/ l ) > ft > diam(r e ). Hence, if -T e nr iDiMil ^ 0, then 
-T e C r xo ^h- From this we deduce that r e C 1^0,5 and T e n A;o,G8-i)/i = for all 
e S I Xo ,s\Ix o ,0h- This implies 

(51) < (/3-l)/i< \x e q -x \ < 5, Ve € I X0>s \I XOtj3h , q = l,...,Q. 

To each quadrature clement _T e with e G I Xo ,s\I Xo ,ph we associate a distinguished 
node q(xo,e) £ {1,...,Q}, radius X(xo,e) £ ((/3 — l)h,S] and subset r% in C T e 
defined by 



mm |x - aj = as (x - x \, 

(52) Q 

X(x ,e) = \x e q{xo ^ e) -x \, r* oM = r e r\B(x ,\(x Q ,e)). 

From (A5) we get Cft < dist(x q ,dr e ) < ft for all q = 1,...,Q, and from this 
we deduce that the area of r% in is bounded from below, namely \r% in | > Cft 2 . 
Combining this result with (A4), we find that 

(53) c\r e \ < \r% o - m \ < \r% y e e i xo ,s\i xo , 0h , e>e^ s . 



18 



O. GONZALEZ AND J. LI 



From flUI} we get, using (A2), (A4) and (A5), and the fact that |Q(^o)| < 1, 
together with ([52]) i and ([53]), 



Q 

e.£l X0 ,s\I X0 ,fih 9=1 

c 



(54) 



< y y — ^^^ e . 



Z^, \rr n _ v e 12— /x 9 



< 



c 

ee/* ,A-r* ./3fc ?(afo,e)l 9 =1 



£ 



— Ixn — X e , s 2 M 



Moreover, by definition of r xoin and the fact that x e q ^ Xo e ) £ dB(xo, X(xq, e)), we 
get 



(55) , r < r, Vy 6 rt 



Combining ([55]) with ([54]) . and using the fact that F Xq in C r e , we find 



|S 2 (s )|< £ I | cL4 y , 



( 56 ) " £ /„. dA »' 



r* 12/ — 



Jr. n . \y-*o\ 2 " 



where the last inequality follows from a direct estimate of the integral using polar 
coordinates and the fact that the local coordinate maps ip Xo an d 4>xo are Lipschitz. 
Thus the result for § 2 is established. 

For the final term S3 in ([44]) . we notice similar to before that r x0i g-h C r xo j 
and that dist(dr xo ^_ h ,dr xo ^) > h > diam(/ ne ). Hence, if r e n r xo ^-h $i then 
r e C ilo^- From this we deduce that r e n r xo ^-h = for all e ^ / XD ,5- This 
implies 



(57) 0<S-h<\y-x o \ Vyer e , e # I Xo , s , E > E fiiS . 
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From ((47]) we get, using (A2), (A4) and (A5), and the fact that \(q(x )\ < 1 and 
j3h < 5/2, together with ([57]), 

— , Q 

|S 3 (zo)|< E ElC(^)^o,<)^ e |, 

e<ZI*o,s 9=1 

Q 



< v V - w e 



(58) e^o.s 9=1 

Q n 

± E E^ e > 

e-^I^o.s 9=1 

< < c^ -2 , VE > Ep, B , xo e r. 

Thus the result for §3 is established. The desired result stated in part (ii) follows 
from (US), ||5BJ|, dSOl) and (|44|) . 

For part (hi), let (5 e (0, d/3], /3 > 10 and E^a > -Eo be given numbers defined 
as before, so < j3h < 5/2 for all E > Ep t s- Moreover, in view of (A7), we suppose 
that E/3 t s is sufficiently large so that diam(supp(Ch)) < 5/2 for all b — 1, . . . , n and 
E > Ep^s- Since the result in part (iii) trivially holds when £0 = %*> we consider 
without loss of generality any arbitrary points xq,x* € r with < \xq — x*\ < 5. 
Given such points, we consider the following collections of surface elements r e , 
e = 1, . . . , E, defined in the same way as before: 

(59) Ix ,/3h C I Xo ,25 C I XOj d and I x ,,ph C I x «,35 C I x ,,d- 

By design, each of the above sets is non-empty, as well as each of Ix ,d\Ix ,2S, 
I Xo ,2s\I Xo ,ph, Ix,,d\Ix»,38 and I x ,,3s\Ix,,i3h for all E > E 0<s . Moreover, we have 
the inclusion I X0} 2S C I X »,3S for all E > Ep^$. Analogous to before, we consider the 
decomposition 

n 

(60) 2J F b(xo) - F b (x*) = Ji(x ,x*) + J 2 (a;o,a;*) + ^{xo, x*), 
6=1 

where Fb(x) = (b(x)H(x,Xb)Wb and 

Q 

(61) J 1 {x ,x,) = E E^^o)-^^*)' 

e£4 ,2i 9=1 

Q 

(62) y 2 (aro,a:*)= E E ^C 3 *) ~ ^C*.). 

Q 

(63) 7 3 (*o,x*) = £ E " ^O&O- 

e£Ix ,d 9=1 

To establish the result, we show that, for every e > 0, there exists a (5 € (0,d/3] 
such that |9t| < e, |5F 2 | < e and |£F 3 | < e for all < |x — | < 5 and £7 > Ep t g. 
Results for the finite interval E £ [Eo,Ep g] are a straightforward consequence of 
part (i) and will be omitted for brevity. 
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For the term Ui in (|60|). we use the inclusion I Xo .2S C 1x^,35, and the inclusions 

^ ,/3h C 1*0,25 and I Xm ,ph C /a:»,35, to obtain 

Q Q 

|ji(x ,^)i< X! E 1^0)1+ E EWi 
(64) =( E + E 

+ ( E + E )EiW 

From the definition of F^, and the same arguments in part (ii) that yielded (|50|) 
and (1561) . we find 



(65) |Ji(sso,z„)l < (CW + CS ") + ( Cht * + C6 ") < C6 "> ^ E > E 0,s> 

where < fi < 1 and the last inequality follows from the fact that h < 5/(2/3). 
Thus, for every e > 0, we can choose 6 £ (0, d/3] sufficiently small to get the uniform 
bound 

(66) \3 r i(xo, x *)\ < £ ! VE > Ep t g, xq^x* £ r, < |xo — x*| < S. 

For the term 3^ in (|60|) . we notice first that -T X0 ,25-h C r xo ^s and moreover 
that dist(0r Bo ,2ff_fc,ar ieo ,; W ) > ft > diam(r e ). Hence, if T e n r xo<2S - h ^ 0, then 
-T e C r xo< 2g. From this we deduce that F e C I^o-d and F e n F Xa ^-h = for all 
e 6 I x „a\Ix ,2S- This implies 

39(5 

(67) ^-<2<5-/ l <|y-a; |<d, Vy £ f e , e e I^a, £ > % 
Thus, by a slight generalization of Lemma |4.2f iii). we have 

(68) \H(x*,y) - H(x ,y)\ < ~ 310 ' Vy e T e , e e I« ,dV*o,M, E > Ep, s . 

\y - zo| d p 

Furthermore, since diam(supp (£;,)) < 8/2, we have 

(69) Q(x ) = 1, C g e (z*) = 1, er e ,ee 4 0) A^o,25, £ > ^,5- 

Combining (|69|l and ([68} with (|62[) . and using the fact that < \xq — x*| < 5, we 
get 

— , Q 

|^ 2 (x ,x,)|< E ElQ^)^^'^)^ -Q^*)^*,^)^!, 



(70) < J2 J2\ n ( x °> x $- n ( x *> x W w z> 

ee/x ,d\-fx ,2« 9 =1 



< 



cs 



y v — ^-wi. 

L ^ Z ^ l^e ^ 



ee/x ,d\-fx ,2« 1- 



=1 I""? 



x" - x 
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Using the same arguments as in (f5"4")) -~(f56" )) with analogous quantities q{xo, e), X(xq, e) 
and F Xq i„, we obtain 

f C 8 

|* a (*o,*.)l< £ j r .]^r^ dA v> 

eeIx ,d\Ia, ,2S 

f OS 
(71) < / i — dA 



< 



CS + CS^, 0<^<1 
CS + CSlnS, n=l 



VE > E, 



where the last inequality follows from a direct estimate of the integral using polar 
coordinates and the fact that the local coordinate maps tp xo and ip~* are Lipschitz. 
Hence, for every e > 0, we can choose 8 € (0, d/3] sufficiently small to get the 
uniform bound 

(72) \? 2 (x ,x*)\ < e, VE>E , 5 , x ,x*er, < \x - x*\ < 5. 

For the last term 3^ in (|60)) . we notice similar to before that r xo< d-h C r XQt d 
and that dist(dr x0t d-h, dr x0i d) > h > diam(r e ). Hence, if F e n E Xo ^-h ^ 0, then 
r e C r X(lt d- From this we deduce that r e n r x<3> d-h = $ f° r au e ^ Ix ,d- This 
implies, using the fact that h < (5/(2/3) < d/ (6/3) , 

(73) ^-<d-h<\y-x \, Vy&r e , egI Xo , d , E>E P , 5 . 

Moreover, since < I Xq — X* | < S < d/3, we deduce from the triangle inequality 
that 

39<i 

(74) -^<\y-x*\, Vyer e , e#I ao , d , E>E M . 
Furthermore, since diam(supp(£f,)) < 8/2 < d/6, we have 

(75) C q (x Q ) = l, Q(x*) = l, \/x e q e r e , e # I X0>d , E>E P , S . 

Combining ([75)1 with (|63p . and using the fact that H(x,y) is Lipschitz on the set 
\y — x\ > 39g?/60 by (A2), and the fact that < \x Q — x*\ < 8, we get 

Q 

\2a(x ,x*)\< Y\Q(x )H(x ,x°)W°-Q(x*)H(x*,x°)We\, 
Q 

< ]T *£\H(x ,x*) - H{x*,x e q )\W e q , 

( 7b > e£T« ,„ ?=1 

Q 

^ E J2C\xo-xJW*, 

< c\r\\x -x*\ < cs. 

Thus, for every e > 0, we can choose 8 € (0, d/3] sufficiently small to get the uniform 
bound 

(77) |£F 3 (a;o,af*)| < e, VE>Ep t s, x ,x*er, 0<\x -x*\<8. 

The desired result stated in part (iii) follows from ((77)1. ((72)1 . ((66)) and ((60)1. □ 
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Our next result establishes some important properties of the local polynomial 
R x (z) of degree p > with coefficients {C Xt o, C x , ai , ■ ■ ■ , C Xjai a 2 —a „} defined in 



Lemma 4.4. Let r satisfy (AO) with Lyapunov radius d > 0, H satisfy (A2) with 
exponent < fi < 1, and {r e ,x e q ,W°,Q,Q} satisfy (A4)-(A8). Then: 

(i) For each x G V and n > no, the local polynomial R x : r Xj d ~> M. kxk is 
invariant to the choice of orthonormal basis in T X F ' . 

(ii) R x (z) is continuous in x G r and z G F x ^ for each n > no and d! G (0, d). 

(iii) For every e > there exists N > such that \R x (z)\ < e for all x G r , 
z G r Xj d and n > N . 

Proof. For part (i), let x G -T and n > no be arbitrary, and let £ x (z) and £, x {z) be 
local Cartesian coordinates in r x c i relative to two orthonormal bases in T X T. Then 
there exists an orthogonal matrix L x G M 2x2 such that 

(78) Ix,a( Z ) = L x.aP^, X ,0(z), Vz G r x>d - 

Mor eover, let {C Xj o, C x>ctl , . . . , C x ^ aia2 ... ap } and {C Xy o, C X:Ctl , . . . , C x } de- 

note polynomial coefficients found using the two bases. From P5 ]) -([28 |) and ([75]). 
and the uniqueness of solutions of (|25]1 guaranteed by (A8), and the fact that 
L^ 1 — L^, we deduce 

(79) C Xt o = C Xt o, C Xyai ... aB = L Xi p iai ■ ■ ■ L Xt p tO , s C Xy p 1 ...0 a , s = l,...,p. 

Combining ([75)1 and (|75|) with (f2"4")l . we find that the value of R x (z) is independent 
of the basis. 

For part (ii), let n > no, d! G (0, d), xo G r and zq G r xo .d> be arbitrary, and 
consider any < S < (d — d')/2. Then, by the triangle inequality, we find 

(80) z G r x ,d, Va; G r x0t g, Vz G r z0t g. 

Moreover, let uq G T Xo r be any given unit vector, let v(x) denote the outward unit 
normal to r at x, and consider a basis for T X F defined by 

f \ uo — (iio ■ v(x))v(x) , . . . , . 

81 h (x = f 2 z =Kz xy- 

lite - ("o • K^jj^Wl 

Then, for any a; G -Txo,<5 and 2 G f Zo ,5, the local Cartesian coordinates for z G -T^.d 
are given by 

(82) U a (z) = {z-x)-t a (x). 

From the fact that these coordinates depend continuously on x G F Xo< g and z G 
r zo s, together with Lemmas 14.21 and 14.31 we deduce that the coefficient matrix and 
data vector for the linear system in (|25p depend continuously on x G r xo _$. Fur- 
thermore, by (A8), the unique solution {C x ,o, C XjCtl , . . . , C x , aia2 ... ap } of the system 
also depends continuously on x G r xo g. Combining these results with (|24[) . we find 
that R x {z) depends continuously 01116 r x0t g and z G F ZOj g as required. 

For part (iii), we notice that, by (A8) and the fact that |£ x (z)| — d f° r all ^ G P 
and z G r x> d, it will be sufficient to show that, for every e > 0, there exists an 
N > such that 

(83) |A°|<e, \A^-P«\ < e, Vn>N, x G T, 8 = 1,..., p. 

For brevity, we establish the above bound for A° only, and note that a bound for 
^Pt~P, follows by a similar argument. Using the same notation as in the proof 
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of Lemma f4.3[ let 5 6 (0, d/3], j3 > 10 and Ep t s > Eq be given numbers, where 
Ep t s is sufficiently large such that < f3h < 8/2 and diam(supp((b)) < 5/2 for all 
b = 1, . . . ,n and E > Epj. Moreover, for arbitrary x = xq G r, we consider as 
before the collections of surface elements r e , e — l,...,E, denoted by I XQ ,ph and 
Ix ,s- From (|26p . we consider the decomposition 

~ n 

(84) A° = / H(x ,y)r) Xo (y) dA y - y2( b (x )H(x ,x b )r] Xo (x b )W b , 
Jr 



(85) 
(86) 



b=l 



E 

e=l 

2) 1 (x ) + D 2 (a;o), 



H(x ,y)r] Xo {y) dA y - ^2Q(x )H(x ,x e q )r] Xo {Xq)W^ 

9=1 



where 
(87) 



Di(*o)= E 



© 2 (a;o) 



e£4 



E 



H(x ,y)ri xo (y) dA y - ^2Q(x )H(x Q ,x e q )i lxo (x e q )W^ 

9=1 



H(x ,y)r) Xo (y) dA y - ^2(^(x )H(x ,x e q )r] Xo (x e q )W q 

9=1 



To establish the result, we show that, for every e > 0, we can choose Epj > Eq 
sufficiently large such that |Di| < e and |1>2| < e for all E > Ep : $ and xq e _T. 
For the term D\ in (f86|) . we use the fact that < 1 to obtain 



(89) |Di(z )|< E 



e£4 



r Q 

/ |ff(a:o,tf)| dAy+Y.lQMHixo^^W, 
Jr " 9=1 



Working with the first sum in (|89|) . we have, using (A2) 



(90) E J J H ( x o,y)\ dA y < J 



, a dA y <C5», VE>Ep. s , 



where < [i < 1 and the last inequality follows from a direct estimate of the 
integral using polar coordinates and the fact that the local coordinate maps ip xa 
and are Lipschitz. Working with the second sum in we have, using the 
inclusion I XOt p h C I Xo ,s, 



(91) 



J2^ x °W x °> x V w « 

e6-f^ n ,5 9=1 



( E + E ) E ic^)^^)^ 

ee^o.^h eG-fx o ,5\-fa: o ,0h 9=1 
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By the same arguments used in the proof of Lemma I4.3f ii) that yielded (f50|) and 
(1561). we find 



Q 

( 92 ) E ^\^o)H(x ,x e q )W^\ < (Chf + CS") <C5», V£ > E M , 

where the last inequality follows from the fact that h < 5/ (2/3). Thus, in view of 
(f92|) . (|90| and (f89| . for every e > 0, we can choose a 5 € (0, rf/3] sufficiently small 
to get the uniform bound 

(93) PiMI < c, VE>Ep, s , x er. 

For the term I>2 in (|86p . we notice that r xo ^s-h C -TL^a and moreover that 
dist(ar Ko ,5_ ft , 5P Xo>5 ) > h > diam(r e ). Hence, if ^nP^,^'^ 0, thenT e C r XD!<5 . 
From this we deduce that _T e n r XQt s-h = for all e Zx D ,<5- This implies, using the 
fact that h < 6/(2/3), 

195 

(94) —<s-h<\y-x \, Vyer 6 , e £ 7 Xo , 5 , E > E M . 
Furthermore, since diam(supp (£;,)) < 5/2, we have 

(95) Q(x ) ^ 1. v ^ G T e , e 4 0)5 , E > E PtB . 

Combining (|95l) with (|88p. and using the fact that (x, y) is class C" 1,1 on the set 
\y — x\ > 195/20 by (A2), and that rj x (y) is class C m,x for all x and y, together 
with the quadrature error bound in (A6) for a rule of order £, we get 

r Q 

(96) |35 2 (xo)| < / ^o.yJ^M^-^fltio.O^KJW; 

(97) < C (5 | J r e |/i min(£ ' m+1) < C s |r|/i min (*' m+1 >, VE>Ep, S - 

Here C5 is a Lipschitz constant for the derivatives of the function H(x,y)rj x (y) on 
the set \y-x\ > 19(5/20, and S E (0,d/3\ is fixed such that (JM]) holds. Thus, for 
any given e > 0, we can choose £73, <s sufficiently large (equivalently h sufficiently 
small) to get the uniform bound 

(98) \T> 2 (xo)\ < e, VE>E^ S , x Q e R 

The desired result for A° o follows from (f98 ]> . t[93 |) and (|86 ]> . □ 

The next result shows that the linear operators A and A n (n > no) defined in 
(j9|) and (fT9|) satisfy the collective compactness conditions (C1)-(C3). 

Lemma 4.5. Let r satisfy (AO) with Lyapunov radius d > 0, G satisfy (Al) } H 
satisfy (A2) with exponent < fi < 1, and {r e , x\, W|, Q, Q} satisfy (A4)-(A8). 
Then: 

(i) A is a compact operator on C°(r, M. k ). 

(ii) A n is a finite-rank (hence compact) operator on C°(r,M. k ) for each n > uq. 

(iii) A n satisfies (C1)-(C3) and hence is collectively compact on C°(r, R k ). 

Proof. The result for A = S + 5£ in part (i) is classic and follows from the continuity 
of G and the properties of H established in Lemma 14.21 see for example |30[ [38] . 
The result for A n = S n + J^n in part (ii) is analogous and relies on Lemma | 
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For brevity, we illustrate only (iii). Moreover, we show the result only for "K n . The 
result for Q n is similar and more straightforward due to the continuity of G. 

To establish (CI), let v € C°(r,R k ) be arbitrary. We seek to show that 
\(3-Cv)(xq) — (3-C n v)(xo)\ — > as n — > oo uniformly in xq G F, where 



(99) 



~ n 

= / H(x 0l y)v(y) dA y -~y2( b (x )H(x a ,x b )v(x b )W b 
Jr 

n 

- E Cb(xo)Rx (Xb)v(x b )- 



6=1 



6=1 



Using the same notation as in the proof of Lemma l4~3l let S £ (0, d/3], ft > 10 and 
Ep,6 > Eq be given numbers, where -Eg,<5 is sufficiently large such that < f3h < 5/2 
and diam(supp(£i,)) < <5/2 for all 6 = 1, . . . , n and -E > £73,5. Moreover, for any 
given £0 S F, consider as before the collections of surface elements -T e , e = 1, . . . , E, 
denoted by I Xo ,ffh and Ix ,6, and consider the decomposition 



(100) 
where 

(101) 3Ci(ar ) 

(102) X 2 (x ) 



{Wv)(x ) - {M n v)(x ) = + X 2 (x ) - X 3 {x ), 



eel. 



E 
E 



Q 



H(x ,y)v(y) dA y - £ Q (x )H(x , x e q )v(x q )W* 



9=1 



H(x ,y)v(y) dA y - ^2Q(x )H(x , x q )v(x e q )W* 



9=1 



(103) X 3 (x ) = >^ Cb(xo)Rx (x b )v(x b ). 



6=1 



To establish the result, we show that, for every e > 0, we can choose Ep t $ > Eq 
sufficiently large such that |3Ci| < e, \%2\ < e and \X 3 \ < e for all E > Ep^ and 

x e r. 

For the term %i in (|100p . the same arguments used to establish (|89 |) -(|92 |) in the 
proof of Lemma 14.41 can be applied to obtain 



(104) 



|3Ci(*o)|< E IM 



\H(x ,y)\ dA y + J2 \C q {x a )H{x Q ,x e q )W e q 

9=1 



<C\\v\\5^ VE>E 0! s- 



Thus, for every e > 0, we can choose a S € (0,d/3] sufficiently small to get the 
uniform bound 



(105) 



|Ki(zo)| < e, VE>Ep, 5 , x e E 
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For the term %2 in (|100|) . the same arguments used to establish ([94 |) -(j97 |) in the 
proof of Lemma 14.41 can be applied to obtain 

Q 



(106) \X 2 (x )\ < £ 



H(x ,y)v(y) dA y - £ H(x ,x°)v(x e q )W< 



(107) < J2 \r e \T(e,f X0 ,h) < \R\ max r(e,f xo ,h), WE > E 0>s . 

Here f Xo {y) = H(xo,y)v(y) is the function being integrated, T(e,f Xo ,h) is the 
quadrature truncation error for this function as given in (A6), e I Xo ,s implies \y — 
xq\ > 19(5/20, and S G (0,d/3] is fixed such that (jl05[) holds. Since H(xo,y) is class 
C" 1 ' 1 in xq and y by (A2) on the set \y — Xo\ > 195/20, it follows that the continuity 
properties of f x „ (y) in y are uniform in xq. Hence, by (A6), max e ^/ x s r(e, f Xo , h) — > 
as h — > uniformly in xq. From this we deduce that, for any given e > 0, we can 
choose Ep : $ sufficiently large (equivalently h sufficiently small) to get the uniform 
bound 

(108) 3C 2 (^o)| < e, VE>E 0! g, x Q e R 

For the term 3C 3 in (|100J . we consider the index set J Xo = {b \ Cb{xo) > 0}. By 
(A7), we have diam(supp(£{,)) < Ch, and by (A5), we have \x a — Xb\ > Ch for any 
a 7^ b. From this we deduce that the number of elements in J Xo , denoted by |J Xo |) 
must be bounded uniformly in xq, namely 

(109) \J X „\<C, VE>Ep, 5 , XQ er. 

From (|103[) we get, using the definition of J Xo and the fact that < £& < 1, 

n 

(110) |3C 3 (z )| ^ J2 \Cb(xo)Rx (x b )v(x b )\ <\\v\\J2 \ R *o(Xb)\. 

6=1 6eJx 

Using (|110[) and ()109|1 . together with Lemma f4.4f iii). we find that, for any given 
e > 0, we can choose Ep_s sufficiently large to get the uniform bound 

(HI) \OC 3 {x )\<e, VE>E fi>Sl x e R 

The desired result follows from ([TTT]) . (fTUSj) . (ITU5|) and (fTUU|) . 

To establish (C2), we seek to show that ( l K n v)(xo) is uniformly bounded for all 
n > n 0l x Q € R and v £ C°{RM. k ) with \\v\\ < 1. By definition of ~K n , we have 

n 

(112) (3i n v)(x ) = ^ \Cb{xo)H(xo,x b )v(x b )Wb + ( b {x Q )R Xo (x b )v(x b ) , 

6=1 

and from the definition of J Xo given above, and the fact that < Q b < 1, we get 

n 

(113) \(X n v)(x )\<\\v\\J2\( b (x )H(xo,x b )W b \ + \\v\\ £ \R Xo (x b )\- 

6=1 beJx 

From this we deduce, using Lemma l4.3f ii) and l4.4f iii). the uniform bound 

(114) \(H n v)(x )\ < C, Vn>n , x G R, veC°(R,R k ), \\v\\<l. 

To establish (C3), we seek to show that (^K n v)(xo) is uniformly equicontinuous 
for all n > no, x Q e T and v e C°(r,M. k ) with \ \v\\ < 1. For arbitrary x ,x* e R, 
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we have, by definition of 

(3< n v)(x*) - (JC n v)(x ) 

n n 

= y^(b(z*)H(x*,x b )v(x b )W b -y^( b (x )H(x ,x b )v(x b )W b 

(115) ti t[ 

n n 

+ ^2Cb{x*)Rx, {x b )v(x b ) - ( b (x )R Xo (x b )v(x b ). 

6=1 6=1 

Similar to before, we consider the decomposition 

(116) (Hnv)(x*) ~ (Wn«)(aro) = Q(x ,x*) + K(ar.) - K(ar ), 
where 

(117) Q(x ,a:*) = [0,(x»)i2"(a;», z&)W6 - C6(xo)#(x , x b )W b ] v(x b ), 

6=1 

n n 

(118) R(x*) = *^2( b (x*)R Xt ,(x b )v(x b ), R(xq) = y^Xb( x o)Rx (xb)v(xb)- 

6=1 6=1 

To establish the result, we show that, for every e > 0, we can choose a S G (0, d/3] 
such that |Q| < e and < e for all x ,x* G T with |x - x*| < <5, u G C°(r,M fe ) 
with 1 1 < 1. and n > hq. Notice that this result immediately follows from Lemma 
I4.3IT iii^) . Lemma l4~4T iii) and (|109[) . after increasing the value of no if necessary, since 



(119) \Q(x ,x*)\ < Mj^Mx^Hix^x^Wb - ( b (xo)H(x ,x b )W b \, 

6=1 

(120) \X(x*)\<\\v\\J2 IKxA^l \0l(x )\ < \\v\\ \Rxo(*b)\. 

6eJx* beJx 

□ 

The next result will be crucial in establishing the linear rate of convergence of 
the lowest-order method with p = 0. We employ the notation from condition (A3). 
Moreover, for any given xo G r and ip G C l,1 {r, R k ), we introduce the function 
Ux .ip ■ r X(ud -> R k defined as 
(121) 

' H(x ,y)[<p(y) - <p(x )] - u»° la *(0,a y ))Dip(x )T Xo (y), y ? x , 

0, y = x - 



U X0 ,Jy) 



Here T Xo (y) G T Xo T is the unit vector defined from the orthogonal projection of the 
chord y — xo onto T Xo T, and £(y) G S is the angular polar coordinate associated 
with y ^ xq. Specifically, we consider the functions T Xo : r XOt d\{xo} — > M 3 and 

-'-xo x o,d 



\y-x \ ~ lx o\y)-> V T x 0, 
0, y = x , 

where vq — v { x o) is the outward unit normal to r at xo- We remark that the result 
in part (ii) (b) below, and its discrete analog in part (ii) (c) , are reminiscent of the 
classic Tricomi condition that arises in the study of singular integral operators [JS] . 
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Lemma 4.6. Let r satisfy (AO) with Lyapunov radius d > and regularity index 
m > 1, H satisfy (A2) and (A3) with exponent fi = 1, and {r e ,x q ,W q } satisfy 
(A4)-(A6). TAen: 

(i) For every xq E T and tp G C 1 ' 1 ^, K fe ) we ftaue U Xo . v G C°' 1 (r x0td ,M k ) 
with Lipschitz constant uniform in xq. Hence \U XOjV (y)\ < C v \y — x$\ for 
all x G r and y G i~x ,d- 

(ii) For every xq G P t/iere ermte a pate/i r* G -T suc/i i/iai 

(a) r XaA / C c r* o c r XOid; 

(b) / r . o U Pf-( 0; £( y ))^(:r )T Ko (y) ^ - for each tp G C^(r,R k ), 

(c) l£°ei; Ei< g <Q^° Iar (0,CK))^K)7;oK)% e | < C v hforalln > 
n and each tp G C 1 ' 1 (_T, I 



In i/ie above, C denotes a fixed constant, C v denotes a constant depending on tp, 
the sum with a prime is over those q such that x q ^ xo, and I* o — {e \ r e G r* Q }. 

Proof. For part (i), let io 6 F and tp G C 1;1 (/ n , M. k ) be arbitrary. By working 
with local Cartesian coordinates in the Lyapunov patch r x0i d and using Taylor's 
Theorem for Lipschitz functions, we find 

(123) p{y) = <p(xo) + D<p(x )[y-x ]+5l Xo (y), Vy e r X0id , 

where 

l^xo (y)l <C v \y-x Q \ 2 , 

Rx £ C 1 ' 1 (-Px ,d)^ )) 1 — ° I G C 0,1 (-Px ,d)^ )• 

Moreover, since 99 is globally Lipschitz on J 1 , the Lipschitz constants for 5l Xo and 
Rx /\y — xq\ are uniform in so £ F, where the latter function is defined to be zero 
when y = xo- 

To establish properties of the functions T Xo and T A , for any given y G r x0i d 
with j/ 7^ xq, wc consider an arbitrary curve 7(7-) G r xa ,d such that 7(0) = y. 
Then, by direct calculation from their definitions, we find |^-T , x (7( r ))lr=o| < 
C|7'(0)|/|y - x \ and \£T*( 7 (t))\ t=q \ < C| 7 '(0)|. From this, we deduce 

\DT X0 (y)\ < j-^— \DT x A o (y)\<C, Vt, G r xo4 , y + x Q , 

(125) \V x o\ 

^gC - 1 ^,* 3 ), 

where the last result above follows from the fact that T^(y) has a uniformly 
bounded surface derivative for almost every y G r XQ ^. Notice that the Lipschitz 
constant for T A , equivalently the bound on its surface derivative, is uniform in 
x G r. Moreover, by definition of T XQ1 and the Lipschitz property of T A , and the 
fact that T*(x ) = 0, we have the bounds 

(126) \T X0 (y)\ = l, \T£(y)\ < C\y - x \, My G r xo4 , y ± x . 
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Proceeding with part (i), for any y G r xo ,d with y ^ xq, we use (|124l) . (|123[) and 
([121]) . and (A2) with /x = 1, to write 

(127) = ff^o.l/JbCtf) - V(^o)] -<° lar (0,e(y))^(x o )T ;co (y), 



= u(x ,y) pv(*o)^f + ^Jfj - «£° I "(0,e(i/))Dv J (»o)T !BO (i/). 

In terms of the functions defined in (|122[) and defined in (A3)(iii), this 
becomes 

U XQ , v (y) = u(xo,y)Dip(x )T^ a (y) 

+ u(x ,y)^^-r + (y)D(f(x )T Xo (y). 

Since \u(xo,y)\ < C by (A2), |T Xo (j/)| = 1 by definition, and \Dip(xo)\ < C v by the 
regularity assumption on <p, and \T*(y)\ < C\y - x \ by (|126p . |^ (y)|/|y - x \ < 
C v \y-Xo\ by (pi, and |<(y)| = |u£ (y) - < (*„)! < C|y-x | by (A3) (Hi), we 
obtain the bound 

(129) \U XQ Av)\ <C v \y-x \, Vyer X0)d) y^x . 

In the above, the constant C v depends on <p, but is independent of Xq and y. No- 
tice that the bound also holds when y = xq by definition of U Xo ^. To establish the 
Lipschitz property of this function, for any given y G A; ,d with y ^ xq, we again 
consider an arbitrary curve 7(7") G r xot d such that 7(0) = y. Then, by direct calcu- 
lation using 0, PS]) . dTStb . (A2) and (A3)(iii), we find |£l/xo,*>(7M)|r=ol < 
It'(O) I - From this, we deduce 

\DU Xo ^{y)\ < C v , Vj/£r xo , rf , y 7^. T , 

where, as before, the last result above follows from the fact that U XQtV (y) has a 
uniformly bounded surface derivative for almost every y G r x0l d- Notice that the 
Lipschitz constant for U Xo>v depends on ip, but is uniform in xq G r. Thus the 
result in part (i) is established. 

For parts (ii)(a,b), let x G r be arbitrary and consider a local Cartesian coordi- 
nate map y = jAz (£) e ^x rf with £ G J?^ jCK 2 , and a local polar coordinate map 
y = ^ Q Mp, 0) = VC lar (p, with (p, G ^ C M+ x S and £ - ro (p, £) = p£. 
Introducing the Jacobian factor 

(131) , (f) = ( de t[(5|L)^]) I/3 , 
we have the relations 

(132) dA y = J(f) dA 6 = J (pi) p dp dS£, 

where dA y is an area element in r Xo ^d, dA^ is an area element in fi X[) ,d and dS£ is 
an arclength element on S. By conditions (L2) and (L3) on r, we note that the 
Lyapunov radius d > can be chosen small enough such that 

(133) 1 < J(0 < 2, V£ G n xo , d , x G r. 
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Moreover, if we let L > 1 denote the Lipschitz constant for tp Xo , which is uniform in 
Xo, then a straightforward argument using the definitions of r xo ^g and f2 Xo ,g shows 
that 

(134) D XQiS/L c n xo , s c D SthS , V<5e(0,d], x e r, 

where D XQi g is the closed disc of radius 5 > at the origin. 

Proceeding with (ii)(a,b), we consider the function I : S — >• K defined by 

{j-d/2L nd/2L "I 

y o j{f£)pdp,j j(-p£)pd P \. 

Notice that /(£) is well-defined since D XQid / 2 L C ^x Q ,<i by (|134[) . From (1135[) . we 
deduce J(— £) = /(£), and from (|135|) and (|133|) we deduce, using a change of 
variable and the fact that </(£)/</(£') < 2 for all S ft XQ ,d, 

d 2 ( d l 2L ~ - f d ^ L ~ d 2 

(136) ^2-y Ap0pdp<I(0< J o J(p£)pdp<p, 

where the second integral is also well-defined since D Xod / L C f? Xo ,d by (|134p . For 
any given choice of orthonormal basis in T Xo T, we consider the components £ = 
(£1,62) and £ = (£1,^2), and the decomposition S — S + U 5_ C K 2 , where S+ and 
5 1 - are disjoint subsets defined as 

S+ = {£eS|£ 2 >0or£ = (l,0)}, 
[ ' S_ = {£eS'|f 2 <0or£ = (-l ) 0)}. 

Moreover, we consider functions i?+ : 5+ — >• M. and i?_ : 5_ — 5- K defined implicitly 

by 



(138) 



J(p|) p ^ = 7(f), V£eS+, 
J{pt) P dp = I(£), V^eS-, 



and from (1136j) and (| 135[) we note that R+, R- E [d/2L, d/L]. Furthermore, we use 
these functions to define subsets of f2 Xo ,d by 

n* xo ,+ = {£ = n Xo4 1 < P < r+(0, e g 5+}, 
«x ,- = {£ = pfe a , d I < p < R-(£), £ e 

and we note that, since R_ S [d/2L, d/L] and S + U 5_ = 5, we have 

(wo) A^/m c (r2* 0!+ u tf* 0> _) c D XOtd/L , Vx e r. 

Finally, we define a subset of r XQjd by 

(141) C = { y = ^ xo (0 e r XOid 1 £ € r?: 0!+ u n* XOt _}. 

The subset -T* has the properties stated in part (ii)(a,b). Specifically, by (|141|) . 
()140p and f)l 34|) . it satisfies the uniformity condition 

(142) r X() d /2L c r* o c r^.d, Vx e r. 
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Moreover, for any ip € C 1 ' (r, M ), we have 




u x °^(0,Oy))D<p(x )T Xo (y) dA : 



(143) = / / ul° laI (0,OD<p(x )T Xo (OJ(pOpdpdSs 

Js + Jo 

r /■«-(?) ^ ~ 

+ / / Ux f^(0,ODip(x )T Xo (OJ(pOpdpdSz, 
JS- Jo 

where we note that the unit tangent vector T Xo (y) is uniquely determined by the 
angular coordinate € S, and hence we write T Xo (0 in the last two integrals. 
Specifically, for any given orthonormal basis {t Xo ,i > ^0,2} m T Xo r , we have 

(144) T So (e)=|iW+fa*»o,a- 
Combining (fT43"l) with (fT38]) . we get 




< lar (0,e(y))^(^o)T, o (y) <L4„ 



(145) =/ < lar (0,0^(^o)T. (0/(Od5 r 

+ / ^° lar (0, £)D<p(x )T X0 (01(0 dS£. 
Js- 

Since u£° lar (0,£) is an even function on S by (A3)(i,ii), 1(0) is an even function on 
S by (|135p . and T Xo (£) is an odd function on S by (|144p . and since the antipodal 
map £ 1— > — £ is a length- and orientation-preserving map of S" onto S* which maps 
5+ onto we find that the two integrals on the right-hand side of (|145[) cancel, 
which implies 

(146) y^ U P° lar (0,e(y))^(a;o)T :co (y)^ = 0, Vx G T, ip 6 C 1 ' 1 (_T, K fc ). 

Thus the results in parts (ii)(a,b) are established. 

For part (ii) (c) , we use notation similar to that used in the proof of Lemma 14.31 
and let j3 > 10L + 1 and E% > Eq be given numbers, where Ep is sufhciently large 
such that < f3h < d/AL for all E > E^. Moreover, for any given xo £ r, we 
consider the collections of surface elements r e , e — 1, . . . , E, defined by 

(147) r xo = {e 1 r e c r*j and i XOif}h = {e \ r e c r xo , 0h }, 

and note that, by design, the sets I Xo ^h C I x and I* \I Xa .ph are non-empty for all 
> E%. Moreover, for any given ip € C7 1,1 (I^, M fe ), we consider the decomposition 

(148) £ ee/ , o E'i< 9 <q ^o«K))% e = 25i(x ) +5? 2 (a ; o), 
where F Xo (£) = u x f™(0,t)D<p(x )T xo (O and 

(149) St^o) = Ee £ 7 I0l3h E'i< 9 < Q ^o(fK))^, 

(150) 23 2 (x ) = E eeW ^ E? =1 ^„(?K))^- 

To establish the result, we show that the above sums satisfy the bounds |23i | < C v h 
and |2?2| < C v h for all E > E% and xq € -T, where the constant depends 
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on ip, but is independent of xq. Results for the finite interval E £ [Eo,Et] are 
straightforward consequences of the boundedness of F Xo and will be omitted for 
brevity. 

For the term Si in (fT48]) . we use the fact that |i*a: | < C v by (A2) and (A3), 
\I xo ,0h\ < C by 611), and E?=i W, e < C/i 2 by (A5), to deduce 

l S i(-o)l < £ ee ^ h E'r< 9 < Q l^o(C«))l% e 

< C^/i 2 < C v h, WE > E}, x £ r, 

which establishes the result for 23 1. 

For the second term ®2 in (| 148[) . we notice that r xo (p—\)h C r xa ^h and that 
dist{dr xoAP _ 1)h , dr x0:ISh ) > h > diam(r e ). Hence, if F e n r xD!(/3 _ 1)ft ^ 0, then 
-T e C r xo ,ph- From this we deduce that, for all e € I Xo \I Xo ,0h, we must have 
-T e n r xo tp—i)h = 0j or equivalently J7 e (~l fi XOl (0_i)/, = in terms of local Cartesian 
coordinate domains. Since D xo ^_ 1 - )h / L C ^ xo ,(p-i)h by (|134[) . and /3 > 10L + 1 
by assumption, this implies 

(152) 0<10/i<|£|, vee^ e , e e I* \I X0 ,p h , e>e;. 

Moreover, for any two points £, £' G f2 e , with polar coordinates £, — p£, and £' = p'£', 
we have 

■g> ^_ \(p-p'W + p'(e-Q\ 

(153) ^ 

<—, V£,£'Gf2 e , e e i: \i xo ,p h , e>e;, 

where the inequality follows from the fact that diam(J7 e ) < Cdiam(F e ) < Ch by 
the Lipschitz property of the coordinate maps ipx and ip X g- Furthermore, since 
<° laI (piO is class C ' 1 in (p,£) by (A3)(i), and T Xo (£) is class (7°° in f by (IT441) . 
we notice that F Xo (£) = u x ° Iar (0, ^)D<p(a;o)T Xo (£) is class C ' 1 in £ uniformly in cc - 
Hence, for each e G I xa \Ix ,Ph we have, using the fact that > 0, and considering 
each scalar component of F Xo £ R fc separately, 



(154) min F Xo (£) < p „,? " < E«£ (0 ■ 



From this we conclude, by the Intermediate Value Theorem, that for each scalar 
component of F XQ there exists an £ e = p e £ e £ Q e with the property that 

Q Q 

(155) E F *o(£K))^ e = ^o(r)E^ e > Vee 4*o\4 ,^, s>^. 

q=l q=l 

For convenience, we use the same symbol £ e to denote the distinguished point for 
each component of F XQ . In what follows, we will have need to consider the difference, 
denoted by r*^ A , between the set r* Q and the subset U ee /»^\/ xo l3h r e . From the 
definition of I* Q and I Xo ,/3h, and the fact that diam(P e ) < h, we deduce that 

(156) r x *> A = r; \( (J r) c r^^u^.ft), 

where N(dr* o ,h) denotes the set of all y £ r* a such that dist(y,dr* o ) < h, that 
is, N(dr* o , h) is the closed neighborhood of dr* Q in r* of size h. We note that 
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r Xo .j3h and N(dr* o , h) are disjoint for E% sufficiently large. Moreover, from (|138[) . 
(|139[) and (jl41[) . we deduce that the set dT* o is a curve of class C° :1 . 

Proceeding with the term 232 in (| 148|1 . we have, using (|155j) and (|150|) . and 
considering each scalar component separately, 



® 2 (* )= E F xo(r)E% e 

= E F *o(r)in + E ^o(n(E^ e -i re 

Using the fact that f r , F Xo (£(y)) dA y = by pl6| . and the definition of r£' A in 

^0 

(j!56|) . we get 

23 2 (x )= e ^(om + E F * r)(E% e -i re 



ai 



( 158 > = E (/ e F* (n-F x My))dA y ) 

e ^ I S \ I x ,Ph 

Q , 

+ E F * (h(J2w°-\n) - l A F xo g(y))dA y 

eei; V* ,i>h 9=1 
= 2 2l (*o) + £ 22 (x ) - ^23(^0), 

where !B2i(xo), 2322(^o) and 23 2 3(a;o) denote the three terms in the middle equation 
of (USB). 

For the term CB21 in (|158|) . we notice that, by the Lipschitz property of F Xo (£) 
and (|153[) . we have 

l^o(n--^o(f(j/))l <C v \t e -tty)\ 

(159) <r C ^ <r C x> h w c r* \ r f > p* 

— — 1 r, ve € I Xo \Ix ,ph, E>E p , 

p{y) \y-x \ 

where the last inequality follows from the Lipschitz property of the coordinate map 
ip Xo and the fact that \y - x \ = \ip Xo (C(v)) ~ ^ (°)l < C\£(y)\ and p(y) = |£(y)|. 
Hence, for the term 2321 we obtain the uniform bound 



|S 2 i(xo)|< E / \F* (n ~ F X MV))\ dA 
e e / i \ / x ,3h 



(160) 

< C v h I - — l — dA y < C v h, VF > E* x e r, 
Jrz \y-%o\ p 

where the last inequality follows from a direct estimate of the integral using polar 
coordinates, which shows that the integral is bounded, and the constant C v depends 
on ip, but is independent of xq. Thus the result for 23 2 i is established. 
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For the term 2322 in (|158[) . we first use the quadrature error bound in (A6) for a 
rule of order £ > 1 applied to a constant function to obtain 

Q Q , 

(161) xx-in = \J2 W «- / dA v ^ c \ re \ h > Ve e % v*o,eh, e > e}. 



Using the above, together with the fact that |-Fz (£)l < ky (A2), (A3) and the 
definition of F Xo , we then obtain the uniform bound 



E 

9=1 



w e q - \r e 



i23 2 2 (^o)i< E 

(162) ee/; \/„ DlPh 



< 



where the constant C v depends on <p, but is independent of xq. Thus the result for 
S22 is established. 

For the term 2323 in (|158p . we again use the fact that |-Fa:o(£)l — @<p as above, 
and use the inclusion result for r*^ A in (|156JI . to get 



|B 23 (aro)| < 



(163) 



/• ; • 



\Fx (Z(v))\ dA y 



< a 



cL4„ + C„ 



: ,/3h 



JV(ar* o> h) 



where the last inequality follows from the definition of the patch r xo ^h and the set 
N(dr* g , h), and the constant C v depends on ip, but is independent of xq. Thus the 
result for 23 23 is established. The result for 23 2 follows from (|TB3"]) . (fTMj) and (|M)1) . 
and the result stated in part (ii) (c) follows from these, together with (|158p . (|151|) 
and (fT48j) . □ 

4.3. Proof of main result. Here we combine the collective compactness result in 
Theorem 14.11 and the results in Lemmas 14.21 - 14.61 to establish our main result in 
Theorem 13.11 We consider the locally-corrected Nystrom method defined in (fTSJ— 
(|28p with a quadrature rule of arbitrary order £ > 1, a local polynomial correction 
of arbitrary degree p > 0, and a surface with regularity index m > 0. 

For part (i), consider any I > 1, p > and m > 0. Then, by Lemma [4.5[ the 
operators .A and A n (n > uq) satisfy the collective compactness conditions (Cl)- 
(C3). Hence, by Theorem 14. 11 there exist constants C v > and N v > n such that 
ap n — A n tp n — f is uniquely solvable for tp n , and moreover 



(164) 



\\y n -v\\<C v \\A n y-Aip\l Vn>N v 



By (CI) and the fact that ip € C°(r,R fe ), we have (A n f)(x) 
in x S which in view of (|164|) implies 



(A(^)(x) uniformly 



(165) 



IK 



as n-^oo, W > l,p > 0, m > 0. 



Thus the result in part (i) of the main theorem is established. 

For convenience, we next consider part (in), and will consider part (ii) afterwards 
since it requires a special treatment. Accordingly, consider any £ > 1, p > 1 and 



A CONVERGENCE THEOREM FOR A CLASS OF NYSTROM METHODS 



35 



m > 1, and assume ip € C m ' 1 (r,M k ). In view of (|164[) . we seek a bound for 
— Acp\ |. By definition of A n and A, we have 

(166) \\A n <p-A<p\\ < \\9 n tp-9<p\\ + \\Hntp-Mp\\, 



where S and 3~C are the operators defined in (fTOj) . and 9n and 3~C n are the operators 
defined in (|18j) : we continue to use subscripts n — > oo in place of h n — > 0. 

For the first term in (|166j) . we let xq G r be arbitrary and consider the func- 
tion (G Xo ip)(y) — G(xo,y)<p(y). Since <p £ C m ' 1 (P,R fe ) by assumption, and G £ 
C"M(r x r,R kxk ) by (Al), it follows that G Xo ip £ C*™' 1 ^, R fe ), where the Lips- 
chitz constant for each derivative is uniform in xq. In view of the quadrature error 
bound in (A6) for a rule of order £, we get, by definition of S and 9«, 

(167) 

E Q 

\(9 n p)(x ) ~ (5p)(x Q )\ = \YY G ^ X >( X ^ W 9 ~ / G(x ,yMy)dA l 

e=l q=l J r 

E Q 

< Yl I Y(G Xol p)K)K ~ I (G X0 <p)(y) dA 



<C v h min ^ m+1 \ Vn>N v , x Q er. 

In the above, the constant C v depends on <p, but is independent of xq. From this 
we deduce the bound 

(168) \\S n <p -9<p\\< C ip h m ^ m+1 \ Vn > N v . 

To bound the second term in (jl66[) , we let Ti o (p : r xo<d R k be the local Taylor 
polynomial of degree j for ip at x S r. Just as in (|24|) . the local polynomial 
(T^. Q p)(z) is defined using local Cartesian coordinates £r (z) in r xo ^- Specifically, 
we have 

(169) (^roV'X 2 ') = G Xo ,0 + Cx ,QiCa:o,Qi ( z ) + Cxo,atia 2 £xo,ati \ z )£xo,ct2 i z ) 

+ ' ' • + G Xo ,uia2---aj(,x ,ai (^)Cxoi«2 ' ' ' £,Xo,otj i z ) ; 

where 
(170) 

r- _ d{ipotp xo ) _ 1 y(yofe ) 

W ,0 — VA^O,), Wo, on — ) ■•• ^xo,ai-aj — ., _,Q£ aj ■ 

In the above, all derivatives are evaluated at the surface point Xq, or equivalently 
£,x {xo) = in local Cartesian coordinates. Since ip £ C m,1 (r, R fc ), it follows 
that T^ip is well-defined for any < j < m. Moreover, in view of the Lipschitz 
properties of the coordinates tj,x (z) for a surface r of class G m+1,1 , it follows 
that T^ o (p £ C m+1,1 (r! EOi d, R k ), where the Lipschitz constant for each derivative is 
uniform in xo. Furthermore, by definition of the locally-corrected operator 3~C n , the 
relation in (|23|) holds, which implies 



(171) {% n r lxo Ti v){x Q ) = {Xr lxo Ti y){x Q ), V0<j<p, x Q £ T. 

Here r] Xo £ C m,1 (r, [0, 1]) is any given cutoff function as described earlier. Notice 
that i] Xo Ti o tp is defined at all points of the surface and rj Xo T^. o p £ C" l,1 (r,]R fe ), 
where the Lipschitz constant for each derivative is uniform in xq- 
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In what follows, we will have need to consider the function ip Xg £ C m,1 (r, R fc ) 
defined by 

(!72) <Px (v) = [<p- Vx T£ <p]{y), j = min(p,m). 

From the definition of the local Taylor polynomial T£ ip, and the fact that r/ Xo is 
identically equal to one in a fixed neighborhood of xq, we deduce that 

(1?3) K(y)l<^|y-a:o| i+1 , \D s ^ o {y)\ <C^\y ~ x^ 1 ^ , 

o<s<j, x ,yer, 

where D s ip x denotes the surface derivative of order s. For any given xq, we will 
also have need to consider the function (H Xo (p^ )(y) = H(xo,y)ip XQ (y). Using (|1T3[) . 
together with (A2) and (A3), we find 

(1?4) \{H X0 ^ ){y)\ <C v \y- ZoF, \D s {H Xa ^ a ){y)\ < C v \y - x»r\ 

o<s<j, x Q ,yer. 

Here and above C v is a constant depending on ip, but is independent of xq and y. 
From (|174p with s = j , we find that the surface derivative of H Xo ip Xg of order j is 
uniformly bounded, namely \D^ {H Xo ip XQ ){y)\ < C v for all x ,y £ r. Thus it follows 
that H XQ if XQ £ &~ 1,1 {r, M k ), where j = min(p, m) and the Lipschitz constant for 
each derivative is uniform in xq. This last regularity result provides a hint of the 
delicate nature of the case p = that will be considered later. 

Proceeding with the second term in (|166j) . we consider an arbitrary xq £ T and 
use ([TTTjt and (fT72]) to write 

(175) (K n ip - X<p)(x ) = (M n <p* - JC^)( Io ). 

Using the definitions of 3~C and "K n in (|175[) . and the fact that C, b + C, b = 1, we get 
(OinV ~ "K<p){x Q ) 

n n 

= (b(xo)H(x , x b )<p£ (x b )Wb + ^ (b(xo)Rx (xb)(p Xo (xb) 

H(xQ,y)<Px (y) dA y , 

= y2,H(x ,x b )ipx Q (xb)W b - / H(x ,y)(p Xo (y) dA y 

n n 

+ ^2(b{x )Rx {xb)v Xo (x b ) - ^2( b {xo)H(x , x b )(p£ {x b )W b . 

6=1 6=1 

For convenience, we decompose the above into three parts, namely 



6=1 6=1 

— / " I -I'D - <I!--' A 

(176) 



(177) (Din<p - Q0p)(xo) = Mi(x ) + M 2 {x ) - M 3 (x ), 
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where 



n „ 

(178) Mi(ar ) = Y / H(x ,x b )<p*(x b )W b - / H(x ,y)<p*(y) dA y , 

n 

(179) M 2 (x ) = Y / (b(xo)R Xo (xb)<fi(xb), 



6=1 

n 

(180) M 3 (x ) =Y,Cb(xo)H(x Q ,x b )<p^ (x b )W b . 

6=1 

To establish the required result, we proceed to bound each of Mi, M2 and M3 
uniformly in Xq. 

For the term Mi in (|T77|) . we note that H Xo ip^ Q <E C^ 1 ' 1 (r, R fc ) by the discussion 
following (|174|) , In view of the quadrature error bound in (A6) for a rule of order 
£, we get 

E Q 

l M i(*o)l = I EE^^K 1 ^ ~ / (^oO(y) dA y 

e=l q=\ ^ r 

<E E(^o)kk- / (h X0 v4)(i/)^ 

e=l q=l 

where the last equality follows from the fact that j = min(p, m). In the above, the 
constant C v depends on tp, but is independent of xq. Thus a uniform bound for 
Mi is established. 

For the term M2 in (|177|) . we consider the index set J XQ = {b | ( b (xo) > 0}, and 
we note by the same arguments that led to (|109[) . that the number of elements in 
this set is uniformly bounded, namely \J Xo \ C for all n > N v and x$ G r. From 
(I179P we get, using the definition of J Xo and the fact that < ( b < 1, 

n 

(182) |M 2 (z )| < 10.(^0)^0(^X^6) I < E l^^llKWI- 

6=1 6GJx 

Since |i4 M| < C by Lemma Eyflni), |^ (^)| < C^N - :roP +1 by and 
£6 — a^o C^ 1 for all 6 £ J Xo by definition of J Xo and (A7), we get 

(183) |M 2 (x ) <C^ +1 = C^h m ^ p+1 ' m+1 \ Vn>N v , x Q e T, 

where the last equality follows from the fact that j — min(p, m). As before, the 
constant C v depends on (p, but is independent of xq. Thus a uniform bound for 
M2 is established. 

For the term M3 in (|177p . we again use the definition of J Xo and the fact that 
< ( b < 1, together with (TTSD]), to write 

n 

(184) |M 3 (x )| < ^2 \tb(x )H(x , x h )<p%{x h )W h \ < £ |(ff l0 <)(^)|ff t . 

6=1 &eJx 

Since \(H XQ ip Xo )(x b )\ < C v \x b - x a \ j by (j!74[) . and |x h - x | < C7i for all b e J Xo 
by definition of J Xo and (A7), and W b < Ch 2 by (A5), we get 

(185) |M 3 (x )| < C v h j+2 = C v h m ^ p+2 < m+2 \ Vn>N Vl x e T, 
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where the last equality follows from the fact that j = min(p, m). Again, the constant 
C v depends on ip, but is independent of Xo- Thus a uniform bound for M3 is 
established. Combining ([155]) . ([I8"5|) . ([151]) and (ITTTj) . we obtain a bound for the 
second term in (|166[) . namely 

(186) WK n y - JOp\\ < C v h m ^ Lp -' m \ Wn > N v . 
Combining (fT86]) . (fT68]) . (fTBe]) and (HMJ), we get 

(187) ||p„ -tp\\ < C ip h min(e ^ m) as n->oo, W > l,p > 1, m > 1. 

Thus the result stated in part (iii) of the main theorem is established. 

For part (ii) of the main theorem, we assume p = and consider any £ > 1 and 
m > 1. In view of (I164p . we again seek a bound for ||./l n y> — Aip\\ and consider the 
decomposition in (|166l) . The term 1 1 S^tp — S^l can be bound exactly as in (|168l) . 
Hence we focus on the term ||3~C„<p — 2£<£>||, which requires a different treatment 
than before due to the limited regularity of H Xo ip^ in the case p = 0. To this end, 
we consider an arbitrary xq £ r and let r* o and I* Q be as in Lemma l4.6f ii). Then, 
using (|175j) and the definitions of 5f and 5f„, we get 

n n 

= ^(b{xo)H(xo,x b )(Px Q (x b )W b + ^( b (x )R Xa (x b )(Px ( x b) 

b=l b=l 



Q 

(188) = 



E 



H(x ,y)tp XQ (y) dA y , 
J2Q(x )H(x ,x e q )^ o (x e q )We - f H(x 0l y)<pi(y) dA y 

.9=1 JrC 

Q . 

^(;(i„)i/(i 0lI ;)<(i;)f;- / 

n 

+ J2(b(x )Rx (xb)tPz (x b )- 
b=l 



For the first term in the first bracketed expression in (|188j) . we observe that 
Q(xq)H (xq, Xq)tfx (Xq) = if Xg = Xq by j 174ft and (A7). Thus, using the notation 
from Lemma 14.61 and the fact that £| + £| = 1 , we have 



J2C q (xo)H(xo,x^ (x* q )W* 

9=1 



(189) =E 1<9 < Q ^o)#(*o,^>^ ; 



E 1 < fl <o ff (»0 > ^)^xo(^)^-E 1 < fl < Q S(»o)H(x 0s X»)<(4)W» 
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Using the definition of U Xo ,<p given in f|121|) . and the definition of ip Xo given in (j 1 72[) 
with p — 0, we get 

Q 

J2C q (xo)H(x ,x e q )<p£ (x e q )W° 

(190) 



9=1 



E x < g < Q u ^lW e q - Q(x )H(x , 

+ E 1<g < Q < lar (0,fK))^(x o )T xo (^)^. 

Similarly, for the second term in the first bracketed expression in (|188[) . we use 
Lemma POir ii) (b) and (I156p . together with the definitions of U XOtV and ip^ a , to write 

V / H{x ,y)<Px (y) dA y 

(191) = E/ H(x ,y)^ (y)dA y - [ uZf al (0,Z(y))D V (x )T Xo (y) dA y , 
= E / U X0 , v (y)dA y - f AU l°^(0,-£(y))D<p(xo)T xo (y) dA y . 

Substituting ([191]) and (|190| into f[188|) . and using the fact that U Xo>lp (x q ) = if 
= xq by Lemma l4.6f i). we get the decomposition 

(192) (J£„p ~ X<p)(x ) = Wi(xo) + N 2 (a; ) - K 3 (a; ) + K 4 (a; ) + N 5 (z ), 
where 

Q „ 



(193) 



Ni(xo) = 2 

eel* lq=l 



J2u X0 , v (x e q )W* - f U Xo , v {y)dA 



(194) 



(195) 



(196) 



eel' 



u*° iaz (0,tty))D<p(x o )T XQ (y)dAv 



W*o) = E Ek < q C*(*o)#(zo, 



K 4 (x ) = £ 



^C|(%)^o,<X«)^ e - / H(x ,y)^ (y) dA 



.9=1 



(197) 



>r 5 (a;o) = ^2Cb(x ) R Xo (x b )(p£ (x b ). 



b=l 



To establish the required result, we proceed to bound each of Ni, . . . , N5 uniformly 
in Xq. 
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For the term Ni, we note that R e C R* C r Xo ,d for all e G I* g by definition, and 
U Xo .ip G C^'H-^o,^ ) with a Lipschitz constant that is uniform in xq by Lemma 
l4.6f iV In view of the quadrature error bound in (A6) for a rule of order £, we get 

|Ni(*o)| < E E^K)^ 6 - / u *o,M dA v 



(198) eei>\=i Jr 
< C v h min ^ = C^h, Vn>N v , x Q £ R. 

For the term N2, we use Lemma l4~6l ii) (c) and (|163p to conclude directly that 

(199) |Wa(x )| <C v h, Vn>N v , x a G R. 

For the term IM3 , we use the same arguments that led to (|185[) to get 

(200) |W 3 (z )| < C^ min ( p+2 ' m+2 > = C v h 2 , Vn>N v , x G R 
For the term N5, we use the same arguments that led to (|183[) to get 

(201) |N 5 (x )| < C v h mi < p+l > m+ V = C v h, Vn>N v , x G R. 

In the above, the constant C v depends on tp, but is independent of xq. Thus uniform 
bounds for DsTi , N 2 , N 3 and 3Nf 5 are established. 

It remains to bound the term N4 . For every xq G -T, we note from Lemma 
|L6^ii)(a) that the set R* g contains the patch R XOt d/c- Since diam(/ ne ) < h and 
diam(supp(£|)) < Ch, we have, after increasing the size of N v if necessary, 



d d 

(202) 2C ~ C 

Vj/Gr e , e?I* , q=l,...,Q, n>N v , xq G R. 



<-~h<\y-x Q \, Q(xo)=0, C|K> = 1, 



Using the above, and the fact that H(xo, y) is class C™' 1 on the set \y— xq\ > d/(2C) 
by (A2), and ip^ is class C"™' 1 by (|172p , both with Lipschitz constants for each 
derivative uniform in xq, together with the quadrature error bound in (A6) for a 
rule of order £, we get 

I Q r 
\Mxo)\ < E E^»OKK- (H X0 ^ a )(y) dA y 

< C v /i min(£ ' m+1) , Vn > N v , xq G r. 

As before, the constant C^, depends on ip, but is independent of xq. Thus a uniform 
bound for K 4 is established. Combining (|2"D5|| , pPIjl . (POO)) , ([TM| . (fl^Hj) and (fl"9"2"j) . 
we obtain the bound 

(204) Ktp\\ < C v h, Vn>N v . 
Combining pH|) . (TTgg)) . (flgrJj) and (|TM|) , we get 

(205) \\<p n - <p\\ < C v h as n -> 00, W > l,p = 0, m > 1. 
Thus the result stated in part (ii) of the main theorem is established. 
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